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Abstract 



\^ . Applying a microscopically motivated semi-classical Langevin description 

^,.^ I of the linear sigma model we investigate for various different scenarios the 

^ ' stochastic evolution of a disoriented chiral condensate (DCC) in a rapidly 

^J^l ■ expanding system. Some particular emphasize is put on the numerical ve- 

l/^ • alisation of colored noise in order to treat the underlying dissipative and 

^^ . nonmarkovian stochastic equations of motion. A comparison with an ap- 

Qv I proximate markovian (i.e. instantaneous) treatment of dissipation and noise 

On I will be made in order to identify the possible influence of memory effects in 

the evolution of the chiral order parameter. Assuming a standard Rayleigh 
cooling term to simulate a D-dimensional scaling expansion we present the 
Mh. probability distribution in the low momentum pion number stemming from 

r~| ! the relaxing zero mode component of the chiral field. The best DCC signal is 

expected for initial conditions centered around (a) Ri as would be the case 
of effective light 'pions' close to the phase transition. By choosing appropriate 
idealized global parameters for the expansion our findings show that an exper- 
imentally feasible DCC, if it does exist in nature, has to be a rare event with 
some finite probability following a nontrivial and nonpoissonian distribution 
on an event by event basis. DCCs might then be identified experimentally 
by inspecting higher order factorial cumulants 9m, {"m > 3) in the sampled 
distribution. 
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I. INTRODUCTION AND MOTIVATION 

The prime intention for ultrarelativistic heavy ion colhsions is to study the behaviour of 
nuclear or hadronic matter at extreme conditions hke very high temperatures and energy 
densities. One of the major goals, particular at the upcoming RHIC facilities, is to find 
evidence for a new state of deconfined partonic matter, the quark gluon plasma (QGP) 
III]. Besides the confinement-deconfinement transition one also expects a transition of hot 
hadronic matter, where chiral symmetry is being restored. Lattice calculations of quantum 
chromodynamics (QCD) give the belief, that both transitions do occur at the same critical 
temperature T^ at vanishing net-baryon densities |Q. 

The formation of so-called disoriented chiral condensate (DCC) |@ has been considered 
as the maybe most prominent signature for the restoration of chiral symmetry occuring in 
the ongoing evolution of the hot matter from the chirally restored to the chirally broken 
phase. The idea is here that in the course of the evolution of the system from the initially 
(and only transiently existing) unbroken phase with (the order parameter being) (qq) ^ 
to the true groundstate with {qq) ^ the pseudo-scalar condensate (qr'y^q) might assume 
temporarily non-vanishing values. This misaligned condensate has the same quark content 
and quantum numbers as do pions and thus essentially constitutes a classical pion field. The 
subsequent relaxation of this field back to the alignment of the outside vacuum could then 
lead to an excess of low momentum pions in a single direction in isospin space. 

The possible occurrence of a semi-classical and coherent pion field was first raised in a 
work of Anselm ^ but the idea of forming DCC was made widely known due to Bjorken 
^ and Rajagopal and Wilczek Pj. Since then many works appeared on various aspects 
of DCC formation in heavy ion collisions. As the microscopic physics governing the chiral 
phase transition is not known well enough, one typically employs effective field theories like 
the linear cr-model in order to describe this non-equilibrium phenomenon. On the other 
hand the description of quantum field theory out of equilibrium is interesting in its own right 
and thus has given raise to a major attraction for theoretical studies in order to describe 
the evolution of disoriented chiral condensates 0, like e.g. standard Hartree factorization 
or large A'c-expansion methods 0. Usually these considerations assume an initial state 
at high temperature in which chiral symmetry is restored by vanishing collective fields. 
Independent thermal fluctuations in each isospin direction of the 0(4) a-model are present. 
This configuration sits on the top of the barrier of the potential energy at zero temperature, 
so a sudden cooling of the system supposedly brings it into an unstable state. This picture 
is referred to as the quenched situation p. The spontaneous growth and subsequent decay 
of these configurations would give raise to large collective fluctuations in the number of 
produced neutral pions compared to charged pions, and thus could provide a mechanism 
explaining a family of peculiar cosmic ray events, the Centauros 0. A deeper reason for 
these strong fluctuations lies in the fact that all pions constituting the classical and coherent 
field are sitting in the same momentum state and the overall wavefunction can carry no 
isospin [|TU|. 



The proposed quench scenario ||^ assumes that the effective potential governing the 
evolution of the long wavelength modes immediately turns to the classical one at zero tem- 
perature. This is a very drastic assumption as the soft classical modes completely decouple 
from the residual thermal fluctuations at the chiral phase transition temperature. Such an 



idealized scenario of immediate decoupling might hold effectively if the expansion and the 
associated cooling of the fireball occurs sufficiently fast ]ll|. An alternative, the annealing 
scenario [1^, was suggested by Gavin and Miiller. They used the one-loop effective poten- 
tial instead of the classical one including thermal fluctuations. For moderate expansion and 
cooling it was shown that the system can exhibit longer instable periods and thus should lead 
even to a stronger enhancement of the soft pionic fields. On the other side, both scenarios 
assume that the initial fiuctuation of the order parameter at the beginning of the DCC for- 
mation are centered around zero with a sufficiently small width in a rather ad hoc manner. 
Preparing the initial configuration with stronger initial fiuctuations, no DCC formation has 
been observed |12[ . If the soft field remains in thermal contact with the fiuctuations giving 
raise to the one-loop potential, then one also has to allow for appropriate thermal fiuctua- 
tions in the initial conditions |T^,|T5[. The proposed quenched initial conditions within the 
linear sigma model seem statistically unlikely. 

The likeliness of an instability leading potentially to a DCC event during the evolution 
with a continuous contact with the heat bath of thermal pions was investigated by Biro 
and one of us by means of simple Langevin equations |T^. There the average and statistical 
properties of individual solutions were studied with the emphasis on such periods of the time 
evolution when the transverse mass fi± of the pionic modes becomes imaginary and therefore 
an exponential growth of unstable fiuctuations in the collective fields might be expected. It 
was found that for different realistic initial volumes individual events of an ensemble lead to 
sometimes significant growth of fiuctuations [0,1131 • Subsequent investigation by us in fact 



leads to the idea of stochastic formation of DCC for particular special stochastic evolution 
of the order parameter 



This idea is what we want to detail in the present study in more depth. Our main con- 
ception is that the order parameter before and after the onset of the chiral phase transition 
still interacts (dissipatively) with its (nearly) thermal surrounding of thermal (or 'hard') 
pions, which then give raise also to large fiuctuations in the evolution. This one can inter- 
pret as a breakdown of the standard mean-field approximation. Applying a microscopically 
motivated semi-classical Langevin description of the linear sigma model we investigate for 
various different scenarios the stochastic evolution of a single disoriented chiral condensate 
in a rapidly expanding system assuming a D-dimensional scaling expansion |Tl],|l3|,|l^,|l6| . 
Our stochastic description will allow for a systematic recording of the statistically possible 
initial configurations of the order parameter. Furthermore it also describes the nontrivial 
infiuence of dissipation and fiuctuations on the nonequilibrium evolution and the coherent 
amplification on the collective pionic zero mode fields during and after the onset of the phase 
transition. 

It remains to answer the important question of how likely particular nonequilbrium 
evolutions of a statistically generated ensemble will lead to the formation of a 'large' DCC- 
domain, which will depend nontrivially on the initial and subsequent fiuctuations suffered 
by the surrounding in the course of the evolution. For this we calculate the effective pion 
number contained in the pionic collective field emerging by the rolling down of the chiral 
fields to its true vacuum values, which by subsequent emission will be freed as low momentum 
pions. With this number at hand we can make the decision whether accordingly these pions 
can contribute to an experimentally measurable enhancement of low momentum pions and 
thus might provide indeed a signal for the occuring chiral phase transition. As it turns 



out the probability distribution in the pion number contains interesting new information 
for the characteristics of the non-equihbrium evolution stemming from the relaxing zero 
mode component of the chiral field. In the interesting cases the to be expected yield in low 
momentum pions do not follow an usual and simple statistical distribution, but posesses 
large and nontrivial (non-poissonian) fluctuations. The best DCC signal is expected for 
initial conditions centered around (a) ^ as would be the case of effective light 'pions' 
close to the phase transition. By choosing some idealized global parameters for a {D =) 3- 
dimensional, spherical expansion, our findings show that an experimentally feasible DCC, 
if it does exist in nature, has to be a rare event with some finite probability following a 
nontrivial and nonpoissonian distribution on an event by event basis. Comparing with an 
additional incoherent background the fluctuations in the low momentum pion number might 
be revealed in the nonvanishing of higher order factorial cumulants 6m (m > 3). Admittingly, 
we have to say that although we do stress a new physical picture our study has still to be 
seen as a fairly idealized scenario. Nevertheless, we believe that our results are interesting 
in their own right and should serve as a simplified estimate for the nontrivial late dynamics 
encountered in an ultrarelativistic heavy ion collision. 

In the next section H we describe the linear a-model within a Langevin treatment. For 
this we will first summarize the theoretical ideas behind a semi-classical Langevin description 
for the soft (i.e. low momentum) fields in thermal quantum field theory. The hard modes 
are treated as thermal quasiparticles which constitute a surrounding, open heat bath. We 
then discuss the model introduced in [jl6[ in more depth for simulating the evolution of the 



order parameter and the collective zero mode pionic fields. The damping term entering the 
dynamical evolution will be discussed and a systematic recording of the statistically possible 
initial configurations of the order parameter for finite volumes will be given. The final 
equations of motion to be used for the dynamical evolution including a D-dimensional scaling 
expansion for modeling the possible formation of DCC are then stated. As a characteristic 
for describing the 'strength' of a DCC we consider the effective pion number content of 
the evolving domain. Some particular emphasize will then be put first in section |ITl| on 
the numerical realization of colored noise in order to treat the underlying dissipative and 
nonmarkovian stochastic equations of motion. This, to the best of our knowledge, is the 
first numerical treatment of nonmarkovian Langevin equations in thermal quantum field 
theory and might be of importance for other related topics. A comparison with a standard 
markovian (i.e. instantaneous) treatment of dissipation and noise will be made. For the later 
simulations it shows that it is sufficient to consider only the markovian approximation, which 
is numerically much easier to handle. In section |I^ we finally present numerical results of 
the simulation on the evolution of a coherent pionic field. Four different scenarios, annealing 
or quench with initial conditions governed by effective 'light' or physical mass pions, will 
be investigated. We calculate the pion number for a single domain and the distribution 
of the pion number, which are the observables relevant to the experimental detection of 
DCC, and which also will give quantitative predicition on the possibility of forming an 
experimentally accessible DCC. The unusual distribution in the number of low momentum 
pions is further analyzed by means of a cumulant expansion. To be more realistic we also 
take into account an additional incoherent (poissonian) contribution on the production of 
soft pions. Inspecting the resulting distribution for the low momentum pion number we show 
that the high order factorial cumulants can still be large. This provides a new signature 



to identify possible DCC formation. Some conclusions for possible experimental searches 
are drawn. We close our findings with a summary. In appendix ^ we give a microscopic 
derivation of the frequency dependence of the dissipation kernel being employed. Appendix 
describes in more detail our strategy for simulating gaussian nonwhite, colored noise for 
an arbitrary noise kernel. Appendix |C| gives a brief reminder on the cumulant expansion for 
statistical distributions. 



II. LANGEVIN DESCRIPTION OF LINEAR SIGMA MODEL 

In this section we develop in some more detail the Langevin description of the linear a- 
model introduced in [^B|,|TB|. The starting point is the phenomenological Lagrangian which 
is given by 

£ = ^ d,(t>ad''r - \{<l>ar - V^Y + if 00 , (1) 

where (j)a = (cr, tti, 7r2, vts). We employ the standard parameter /^ = 93 MeV for the pion 
decay constant, -m^ = 140 MeV for the vacuum mass of the pion and rrifj ^ 600 MeV for 
the 'mass' of the a-meson. For the three parameters in (|I|) one then finds 

A = ^^^^3=^ ^ 20, v' = f^-'^ = {87MeVy,H = f^ml = (122MeV)l (2) 

The linear cr-model represents an effective chiral theory of the low energy properties of 
QCD. It can be motivated in more theoretical depth from QCD by the modern methods of 
the renormalization group W^. At finite temperature, to leading order in A, the thermal 



fluctuations ((50^(x, t)) of the pions and cr-mesons do generate an effective Hartree type 
dynamical mass giving raise to an effective temperature dependent potential. In the high 



temperature expansion this results in [20 



2 {T:^m^,m^) \ / N + 2\ ^ N=4: A 3 /o\ 

m,n -^ 4l^j^-"2^- ^^^ 

The resulting chiral phase transition is compatible with the expectations of lattice gauge 
QCD calculations [^. There exist also convincing theoretical arguments [^ that the chiral 
phase transition near the critical temperature T^. of QCD with two massless quarks lies in 
the same universality class as an 0(4)-Heisenberg magnet and thus (in this idealized case 
of massless quarks) exhibits a true second order phase transition which can be described 
in the Landau-Ginzburg theory by means of an effective linear a-model. In this sense one 
considers the linear a-model as an appropriate realisation of the chiral behaviour of QCD 
over the whole range in temperature, though the effective parameters near T ^ T^ need not 
really be equivalent to those at T = 0. 

We now address in an intuitive model how one can go beyond the mean field level for the 
semi-classical chiral collective fields. Our main physical conception will be that the order 
parameter and the collective fields before and after the onset of the chiral phase transition 
still interacts (dissipatively) with its (nearly) thermal surrounding of thermal (or 'hard') 
particles. To outline these ideas more conceptually we will first summarize in the following 
subsection the theoretical reasonings behind a semi-classical Langevin description of the 
soft, i.e. low momentum fields. 



A. Equations of motion for long wavelength modes in a heat bath 

One of the recent topics in especially nonabelian massless quantum field theory at finite 
temperature or near thermal equilibrium concerns the evolution and behaviour of the long 
wavelength modes. These modes often lie entirely in the non-perturbative regime. Therefore 
solutions of the classical field equations in Minkowski space have been widely used in recent 
years to describe long-distance properties of quantum fields that require a non-perturbative 
analysis. A justification of the classical treatment of the long-distance dynamics of weakly 
coupled bosonic quantum fields at high temperature is based on the observation that the 
average thermal amplitude of low-momentum modes is large and approaches the classical 
equipartition limit 

IpHo T_ 
m* 

in the case for a sufficiently small generated dynamical mass m* -C T. On the other hand 
the thermodynamics of a classical field is only defined if an ultraviolet cut-off kc is imposed 



nK) = (e^-p/^-l)"^l^°^»l (4) 



on the momentum p such as a finite lattice spacing a. In a recent paper p3[ it was shown, 
at least principally, how to construct an effective semi-classical action for describing not 
only the classical behaviour of the long wavelength modes below some given cutoff kc, but 
taking into account also perturbatively the interaction among the soft and hard modes. The 
resulting effective action 6*08? [soft], which one has to interpret as a stochastic, dissipative 



action p3|,|4|, turns out to be complex, leading to a stochastic equation of motion for the 



soft modes. If the hard modes are already in thermal equilibrium then the evolution of the 
soft modes is described by a set of generalized Langevin equations - the equations of motion 
corresponding to the above complex effective action. 

We briefly sketch the main strategy following [^ by considering a scalar field with 
interaction £int = fr 0'^- The splitting of the Fourier-components, (f){p, t) = (j){p < k^, t) + 
ip{p > kc,t), leads to the following interaction part in the action 

^int[0, ^] = - Ifx (^^^' + 1^ (0V + l^V + ^^')) ■ (5) 

By integrating out the hard modes up to second order in the interaction, one obtains the 
effective action (or influence functional) Sif[4>, (f)'] for the soft modes following the Feynman- 
Vernon approach ||25|. Fig. |l] shows the resulting non- vanishing diagrams contributing to 



SiF- The contributions from diagram (a) and (b) are real and generate the Hartree like 
dynamical mass term. Moreover one notices that Feynman graphs contributing at order 
0{g^) (diagrams (c), (d) and (e)) to the effective action contain imaginary contributions. 
Their real part leads to dissipation (like in linear response theory) whereas the imaginary 
part drives the fluctuations of the hard particles on the soft modes. From the effective action 
semi-classical, stochastic equations of motion result, which have the general shape 

^<P + m^ + C'^^ + E m^^^"''"^^^^-^)^"' = E <P''''^N . (6) 

"J- Ar=i l^^^ ^)- N=l 

Here r2Ar denotes the effective contribution with 2N soft legs, fh? the resummed Hartree- 
Fock self energy (cactus graphs) and ^n are associated noise-variables with a correlation 
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(Cn^n) = Inir2Ar. These generalized Langevin equations (|) are similar in spirit to those 



obtained by Caldeira and Leggett in their discussion of quantum Brownian motion ||2^ . 

For the sake of simplicity we concentrate from now only on the contribution of the 
sunset diagram (c) of Fig. (|l|), i.e. the N=l contribution of (§). Performing a Fourier 
transformation to (H) yields the semi-classical stochastic field equation for a soft mode with 



momenta k ||23|j2 



0(k,t) + (k^ + m*') 0(k,t) + ^J''fh^ e{K - |k - ki - k2|)0(ki,t) X 

0(k2, t)0(k - ki - k2, t)+2 f dt' r(k, t - t') 0(k, t') = e(k, t) . (7) 

J — OD 

r(k, t — t') and ^(k, t) denote the (real valued) dissipation kernel and the noise source, 
respectively, due to the thermal fluctuations of the integrated out hard particles. The 



dissipation kernel is related to the standard imaginary part of the sunset diagram via [24 



Tik,uj) = ^^^, (8) 



which follows by a partial integration of 



/ dt' Sret(k, t - t') 0(k, t') = -2r(k, At = O)0(k, t) + 2 f dt' r(k, t - t') 0(k, t') . (9) 



-oo 

-oo 



(The integration constant represents an additional momentum dependent shift in the dy- 
namically generated mass and will be neglected further on.) The explicit calculation of the 
dissipation kernel is given in the appendix 0. 

Within the present treatment the noise turns out to be gaussian, but colored, character- 



ized by the (ensemble averaged) correlation function [ p3| , p^ 

((^(k, t)^(k', t'))) = {2n)'6%k + k')/(k, t - t') (10) 

or 

((e(k,t)e(-k,t'))) = \//(k, t - t') , (11) 

where the noise correlation strength is related via a generalized fluctuation dissipation rela- 
tion to the dissipation kernel as 

/M=-^^5M±irM'^-2rrM. (12) 

exp(rt,t<;/T) — 1 

In the high-temperature limit 00 <^ T the noise acting on the dynamics of the soft modes 
then fulfills the (entirely) classical relation 

((e(k, t)e(-k, t'))) = 2TVT{k, t-t'). (13) 

The fiuctuation-dissipation-theorem ensures that the soft modes approach thermal equilib- 
rium precisely at the temperature T of the hard modes. 



When the characteristic time scale in the evolution of hard modes in the heat bath is 
short compared to the one of the soft fields and its coupling to the soft fields is sufficiently 



weak, the appropriate ('instantaneous') markovian limit then has a form [E3 



2 / dt'T{k,t-t')<j){k,t') ^ r]<j), (14) 



where rj = T{uJk = \rm?+W) in the linear, harmonic approximation describes the familiar 
on-shell plasmon damping rate (see also appendix ^. In the semi-classical, high temperature 
limit and within the markovian approximation the noise becomes white, i.e. 

((e(k,t)e(-k,t'))) = 2TVv6{t-t'). (15) 



B. Effective description of zero mode in the linear a-model 

In an ultrarelativistic heavy ion collision the idealized onset of a 'quench', as assumed in 
^, is not really be given. Instead, one expects that the most dominant particles to be freed 
after the onset of the transition are the light pions, which represent a thermalized, further 
evolving system. Their occupation in phase space is described via a Bose distribution and 
cannot be correctly taken care of in a purely classical field description. This environmental 
pion gas may then actually expand rapidly enough (in longitudinal and transversal direc- 
tions) to allow for a nonequilibrium rolling down of the chiral order parameter and giving 
potentially raise to the formation of a DCC In any case this gas of 'hard' pions does repre- 
sent a heat bath with which the order parameter and the long-wavelength coherent pionic 
fields do interact. In this sense these collective modes represent an open system, which acts 
dissipatively and fluctuatively with the environment. Referring to the general ideas outlined 
in the previous subsection one thus expects that the (assumed semi-classical) dynamics of 
those modes can be described by means of appropriate Langevin equations. This intuition 
gave the phenomenological basis for the equations of motion used in [^ 



As we will argue in subsection [II D| we expect that for realistic initial (small) volumes 



V{to) the zero mode (k = 0) pionic fields will cover the dominant coherent pion modes to 
be possibly amplified in the course of a sufficient rapid evolution of the system. These three 
modes in fact do represent the pionic portion to the zero mode 

$'^(t) := l.Jd^a:r{x,t) ^ {a,7r\7c',^){t) . (16) 

In the following we want to restrict ourselves to the effective description of this zero mode 
field $", which formally corresponds to the limit k^ ^ in the previous discussion. 

In analogy to (^ we now propose the following effective Langevin equation of motions 
for the zero mode 

% + T[%]+ 121^0 = Uml + ^o, 
$, + r[<i>,]+/il<l>, = ei. (17) 



The temperature dependent one-loop transversal ('pion') and the longitudinal ('a'-meson) 



mass for the respective fluctuations (see also (|)) are given by ||13|-[T5[ 

= A($^ + E^? + ^r'-/.')+-^ (18) 

/^S = /^1 + 2a($^ + E$?)- (19) 

The dissipation functional r[$] as well as the (semi-classical) noise will be treated either in 
the markovian approximation or within the full nonmarkovian expression 

{?]$ (markovian appr.) 

t . , (20) 

2 / dfT{t-t')<^{t') (nonmarkovian) 
— oo 

(ium - , uM(.m - 1 It!"-'':!:^ ''^:'''''""\ . (21) 

[ —I (ti — t2)0ab (nonmarkovian) 

Here T denotes the temperature and V the size of the volume of the considered system. It 
will be the major point of section |T| to simulate nonmarkovian Langevin equations and to 
compare them with the appropriate markovian treatment. 

These coupled Langevin equations (|T^ resemble in its structure a stochastic Ginzburg- 
Landau description of phase transition ||2^, especially for an overdamped situation |]28| 



where the $-term can then be neglected. On the other hand, with A ~ 20 we are obviously 
not in a weak coupling regime, so that the formal apparatus layed down in the previous 



section ^IA\ can only serve as a basic motivation. Semi-classical Langevin equations may 
not hold for a strongly interacting theory as for highly non trivial dispersion relations the 
frequencies of the long wavelength modes are not necessarily much smaller than the tem- 
perature. Still, when the soft modes become tremendously populated one can argue that 
the long wavelength modes being coherently amplifled behave classically Pj. Aside from a 
theoretical justiflcation one can regard the Langevin equation as a practical tool to study 
the effect of thermalization on a subsystem, to sample a large set of possible trajectories 
in the evolution, and to address also the question of all thermodynamically possible initial 
conflgurations in a systematic manner. 

A physically motivated choice for the damping coeflicient and the dissipation kernel F we 
will state immediately below. For the moment we stay to the markovian case and take 77 as 
an appropriate free parameter. The 'Brownian' motion of the soft fleld conflguration leads to 
equipartition of the energy at constant temperature. In Fig. ^we show the effective transver- 
sal masses fi± of the pion modes and /iy of the a mode as a function of the temperature 
obtained by solving eqs. (|1^ at flxed temperature T and sufficiently large volume V. The 



masses shown are thus taken as an ensemble average of the different realizations within the 
Langevin scheme. For larger volumes the fluctuations in the obtained masses are of the order 

9 



\jV and thus small. For the situation that the vacuum pion mass is assumed to be zero (no 
explicit symmetry breaking) one can realize from Fig. ^ the situation for a true second or- 



der phase transition occuring at the transition temperature T = Tc^ V 2/^ — 2m^/\ ^ 125 
MeV. On the other hand for the physical situation of a nonvanishing pion mass of ?7i^ = 140 
MeV the 'phase transition' resembles the form of a smooth crossover. In this case, atT ^ T^ 
the (T-field still posseses a nonvanishing value of (cr(T ^ T^)) ~ fT,/2 ^ a^ac/'^- (In the large 
volume limit one has (a) = (|$|) at fixed temperature, where |$| = V<^^ + tt^ denotes the 
magnitude of the order parameter.) Comparing with results of lattice QCD calculations the 
transition temperature of T^ ~ 125 MeV is considerably smaller than the typical ones of 
Tc ~ 150 — 200 MeV. This one might correct by using instead of (|^) the value obtained by 
the large N-expansion |^ as m^^ = ^T^ resulting in T^ ~ 154 MeV. On a qualitative level 
the present description of the chiral phase transition is compatible with the expectation of 
lattice calculations. However, for the later one finds that the phase transition occurs in 
a much sharper window around the critical temperature T^. Slightly above Tc the order 
parameter (qq) ~ (a) already nearly vanishes; furthermore, sufficiently below Tc, the order 
parameter has merely changed from its vacuum value. This abrupt behaviour around the 
critical temperature is not realized within the present treatment of the linear a-model, which 
obviously shows a much smoother dependence with temperature. A more refined analysis 



within the linear cr-model might account for this behaviour [^ . 

We now turn our attention to specify the dissipation coefficient rj or damping kernel F 
of ( pOD entering the Langevin equations ([T7|). From a physical point of view they should 



incorporate the net effect of the dissipative scattering of the thermal ('hard') pions with the 
collective fields. Its value is thus also of principal interest for DCC formation as a (too) 
large damping of the collective pionic fields would subsequently reduce significantly the 
amplitude of any coherently amplified pionic field [0 and thus might destroy any possible 



DCC. Being consistent within the linear a-model we consider here the 'sunset'-contribution 
(see Fig. |l](c)) as the dominant term for the dissipation, as it incorporates the net effect 
due to scattering of a soft mode on a hard particle into two hard particles and vice versa. 
The on-shell plasmon damping rate can then easily be evaluated in analogy to standard 
0'^-theory to be 



9 ,.T 



2 



'tp 



V=j^^'—fsp{l-e-^) , (22) 

lOTT'^ nip 



where fsp{x) = — /f dt^ (see appendix ^). As emphasized in |2^, the appropriate marko- 
vian approximation in a weakly coupled theory just corresponds to this on-shell approx- 
imation. At first sight, in the present situation of a strongly coupled theory, one might 
think that this 'choice' can only provide a rather crude estimate as the zero mode do not 
evolve on-shell during the (possibly unstable) evolution. Hence the dissipation and noise 
correlation should better be described by nonmarkovian terms including memory effects. 
For this we have to evaluate the complete (off-shell) frequency dependence of the dissipa- 
tion kernel. This calculation we have shifted to appendix ^. As a further assumption we 
now take for the plasmon mass rUp the 'pionic' mass n±{T) for the transversal fluctuations 
depicted in the right upper plot of Fig. |^. This choice should be valid near or above T^ as 
the transversal and longitudinal masses becomes nearly degenerate. The thus resulting dis- 
sipation coefficient r] of ( P^ ) is shown in Fig. ^ as a function of the temperature T (see also 



10 



With this prescription one notes that r] posesses a maximum value of ~ 100 MeV near 
the critical temperature, which will result in relaxation (or equilibration) times of roughly 2 
fm/c (compare also with Fig. |^). For sufficiently smaller temperatures rj decreases fast to a 
neglible small value as the density of the thermal pions as potential scattering centers also 



falls rapidly with decreasing temperature. This behaviour is in line with findings in pO 
where the on-shell damping coefficient has been calculated by means of standard chiral pion 
scattering amplitudes in the vacuum. 

Some critical remarks are in order: (1) It is questionable that above the critical temper- 
ature all contributing degrees of freedom are being considered. Above Tc one expects that 
due to the deconfinement transition occuring at the same critical temperature quarks and 
gluons are freed and thus might have a considerable influence on the damping coefficient of 
the collective, mesonic excitations. (2) The damping coefficient rj introduced in ( p2D should 
be appropriate for temperatures close to T^, where spontaneous symmetry breaking has just 
emerged. On the other hand, for a deeply broken phase (T -C Tc), the n + n ^ % + tt scat- 
tering amplitude will become significantly reduced by the additional t-channel exchange of a 
(T-meson, leading to the well known chiral derivative coupling for lower transferred momenta. 
This additional contribution for the deeply broken phase we have not taken into account and 
we thus overestimate the damping associated with the thermal scattering especially for low 
temperatures (see eg for comparison the damping coefficient given in ||3^). (3) Moreover, for 
temperatures much below T^ the 0(4) transverse and longitudinal mass for the fluctuations 
are not equal anymore. From Fig. |^ one recognizes that below T ^ 100 MeV the longitudi- 
nal mass /i|| exceeds two times the transversal mass ^±, so that the decay of the longitudinal 
mode into two transversal particles becomes possible. In vacuum this just corresponds to 



the decay a -^ vrvr |3^;34] with a width on the order of a few hundred MeV. This would 



give raise to an additional temperature dependent dissipation in longitudinal direction for 
the evolving order parameter and might have also interesting consequences for the DCC 
formation investigated in section ^. Qualitatively one expects that the associated damping 
will then effectively slow down considerably the rolling down in longitudinal (i.e. 'radial') 
direction of the order parameter along the effective potential. We leave an implementation 
of this kind of longitudinal damping for future work. (4) A flnal problem, which we briefly 
mention, concerns the chiral limit rrij^ = 0. Below T^ the pions remain as massless Goldstone 
bosons (see also Fig. ^ and the cr-meson becomes degenerate with the pion at and above the 
critical temperature. Taking the expression (0) one notices that the dissipation coefficient 
7] diverges like l/rup. On the other hand one expects for a true second order phase transition 
a critical slowing down of the excitations near the critical temperature and thus a vanishing 
of the dissipation coefficient ||HT[]. This then implies that a perturbative evaluation is not 
valid but requires a nonpertubative analysis via e.g. renormalization group methods [^ . 



Our discussion should demonstrate that a precise determination of the description of 
the dissipation functional r[$] near the critical temperature is far from being settled. We 
consider our choice as a physical motivation, and which is also numerically tractable. 
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C. Fluctuations of initial conditions at critical temperature 

As a first and straightforward application we address the important question for the 
possible distribution of the order parameter (0) at the critical temperature for a finite 
system with fixed size V. With the noise fiuctuating according to (^) we expect (similarly 
like in Brownian motion) that the chiral fields do fiuctuate thermally around its mean as 
well. Assuming that slightly above the transition temperature the system is near thermal 
equilibrium, generating an ensemble distribution then offers a systematic sampling of all 
possible 'initial' configurations for the later dynamical evolution of the order fields, which 
then lead to a stochastic formation of DCC 

In Fig. ^ we show first the sensitivity of the (ensemble) averaged value of the order 
parameter ((|$|)) = ((vo^T^)) on various sizes V as function of the temperature. As 
expected, finite sizes lead to a positive shift of the order parameter and to a (further) 
rounding of the phase transition. In Fig. |^ the characteristic distributions of the chiral fields 
and their 'velocities' at the critical temperature are depicted. The average width scales like 



l/y/V. Such a behaviour has been reported already within an independent approach in ^E 
One might also employ the quantal version of the noise fiuctuations according to (|12]), which 
in the markovian on-shell treatment one would approximate as 



Trip nip 



yVcoth{^)6{h-h). (23) 

Such a prescription results in even larger fiuctuations. It is also interesting to look at the 
situation in the chiral limit m.^ = 0. The characteristic distribution P{cr) is given in Fig. 
p. In this case the fiuctuations are even larger and scale effectively with 1/V^^*. (One 



can find analytically [|T^ that for this case ((a^)) = 1/2 JTc/ {XV), so that the width in the 



distribution P{a) thus has to scale with V~^^^.) 

In the next subsection we will now turn to the description of the chiral fields for an 
expanding environment leading then to stochastic individual trajectories with considerable 
fiuctuations and thus also for particular events out of an ensemble possibly to experimentally 
accessible DCC candidates. In a sense the 'faith' of all individual trajectories (entering to 
some amount the unstable region /x^ < |]16[) is not really predictable and has to be sampled 



in some quantitative way as within our proposed Langevin picture. We have to admit that 
one can certainly improve in various ways on many aspects in describing phase transitions 
out of equilibrium. Much retains to be learned about how these condensates evolve in 
out-of-equilibrium. Probably the most ambitious description on the quantal evolution of 
the chiral fields in out-of-equilibrium has been developed by Niegawa [^ employing the 



powerful closed-time-path (CTP) real-time Greens function technique. It has to be seen 
whether this formal development can be used for practical simulations concerning DCC 
formation. Using the CTP technique, this approach (as well as earlier developments in the 
same direction [§) is, by construction, an ensemble averaged description |2^, which can 
thus describe within sophisticated methods the dynamical evolution of (ensemble averaged) 
expectation values. Unusual fiuctuations, like e.g. in the pion number, as shown later here, 
can only be accounted for by higher order correlation functions. These are typically not 
considered. Our approach, we believe, states thus a fresh new way in order to account in a 
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simple transparent manner for such unusual strong fluctuations and being far from a simple 
gaussian mean field treatment. 

D. Modelling the evolution of potential DCC 

In the following we will state the final equations of motion for simulating the stochastic 
formation of possible DCC. In the markovian approximation these corresponds to the ones 
proposed originally in |1^. As a later characteristic quantity we will consider the pion 



number of the zero mode contained in the evolving domain, which is assumed to roughly 
correspond to the effective number of soft pions freed from the subsequent decay of the 
pionic fluctuations, i.e. the final decay of the DCC. 

It is instructive to first outline how possible DCCs would be formed in a heavy ion colli- 
sion. This intuitive and idealized physical scenario will give some insight for the choices of 
the value of the free parameters to be specified and will also give a perspective to understand 
the physical matters to be discussed in the following sections. Our picture of a possible DCC 
formation in high energy heavy ion collisions is as follows: 

• In the first stage of the collision (at proper times r = 0.3 . . . 0.5 fm/c in the respective 
subvolume of the system) a parton gas is formed with a temperature T ^ T^ well above the 
chiral restoration point. Chiral symmetry is completely restored in this hot region. 

• In the following (r = 2 . . . 3 fm/c), because of the subsequent collective expansion (longitu- 
dinally or later even transversally) the temperature drops to around the critical one (T k. Tc) 
and some small chirally restored or already slightly disoriented domains of collective pionic 
fields start to form together with a thermalized background of (quasi-)pions and possibly 
other thermal excitations within the respective subsystem. The individual subsystems are 
assumed to evolve independently as they are spatially separated and might be separated in 
rapidity. The possible distribution of the chiral (mesonic) order parameter then depends on 
the size of the volume V{t) of the individual domain, as shown in the previous subsection. 

• At a further time (tq = 3 ... 7 fm/c) the temperature of a (rapidly) expanding domain 
crosses the critical temperature T^, having a certain volume V{tq). At the same time the par- 
tonic gas would undergo the deconfinement /confinement phase transition into the mesonic 
freedoms. The temperature of the surrounding 'heat bath' further decreases as the volume 
increases due to the collective expansion. At this stage chiral symmetry becomes sponta- 
neously broken. The stable point of the order parameter characterizing the broken phase 
moves from (a ~ 0, vf ?« 0) in the symmetric phase towards {a k, f^^jr k, 0) in vacuum. This 
change happens fast if the system expands sufficiently rapidly. A possible (but not neces- 



sary) instability might arise depending on the actual ('initial') values of the order fields [16 
In certain cases, depending crucially on the 'appropriate' initial configuration, the order pa- 
rameter can 'roll down' in a 'disoriented' direction with a fixed orientation in isospin space, 
giving raise to a large coherent collective pion mode. A potential DCC is formed. Possible 
DCC domains differ from each other in the orientation in isospin space, in the size and in 
the pionic content. A large DCC domain denotes here a large pionic content. Intuitively 
the order parameter in such a large DCC domain will go through a trajectory deviating 
strongly from the cr-direction during the roll-down period. In any case a sufficiently fast 
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expansion and cooling is mandatory for the possible formation of larger DCCs. (Because 
of the explicit symmetry breaking term Ha, which, in analogy to a ferromagnet, acts as an 
external and rather strong constant magnetic field, together with the dissipative interaction 
with the heat bath, the order parameter will otherwise align more or less quasi adiabatically 
at its thermally dictated equilibrium value along the a-direction, if the experienced cooling 
is not fast enough.) 

• With the ongoing (radial) expansion (r > 10 fni/c) and due to the explicit symmetry 
breaking the order parameter will oscillate with decreasing amplitude around the stable 
point < a >= /jr along the chiral circle (a^ + vr^ = /^). The expansion will come to a halt 
at some freezeout time, the fireball breaks off. The coherent semi-classical pion state within 
the possible DCC domain decays by the emission of long wavelength pions, with isospin 
distribution characteristic to DCC, and which in number correspond approximately to the 
effective pion number stored originally in the coherent state. If this number of the coherently 
produced low momentum pions is not too small compared with incoherent low momentum 
pions from other, random sources, constituting the 'background', a careful event-by-event 
analysis can provide identification of the DCC formation. 

In the following we want to investigate the evolution of the zero mode chiral fields in 
contact with the heat bath constituted by all the other modes solely in one single domain 
being created out of the initially hot region by means of equations of motion analogous to 
(0). As outlined above, of course, many of such domains might well be created. We assume, 
for simplicity, that these individual domains are independent and do not further interact. 

The (rapid) expansion can be incorporated effectively by means of the boost-invariant 



Bjorken scaling expansion |33] assuming that the order fields $„ = ^ai'T) depend on time 



only implicitly via the proper time variable r = Jt'^ — x^rr, where Xejf '■= z for [D =) 1- 
dimensional longitudinal expansion and Xg// := r for [D =) 3-dimensional radial expansion 
||TT| , p!3|JT^ , p!6|j35| , |37|j38| . In the equations of motion the d'Alembertian is then replaced by 
(9^/(9r^ + {D / T)d / dr , giving raise to an effective Raleigh damping coefficient D /r. This one 
might also interpret as an effective Hubble constant |33] due to the volume dilution 



V D /r \^ 

_ _ _ . _ ^(r) = V(r„) (-) (24) 

for the expanding volume V{t) of the domain. In the quasi- free regime of a freely moving 
bosonic field the amplitude then decreases in (proper) time with ~ r~^/^. 

From (p!?!) we then receive the equations of motion for the zero mode fields in an expand- 
ing environment as 



$0 + - $0 + r[<i>o] + /ii$o = Uml + e, 



r 



0; 



$, + -<i>, + r[<i>,] + /ii$, = e^. (25) 



r 



The dissipation functional r[$] as well as the transversal mass //^ do both depend on the 
temperature T{t). The stochastic noise fields obey (|21|). One therefore also needs to know 
how the local temperature evolves with time. In principle one has to ask for the equation 
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of state of the system and solve for the hydrodynamic equations within the (assumed) D- 
dimensional scahng expansion. For the ideal case of a massless gas (which is not a too bad 
approximation for pions) an isentropic expansion results in 



Tir) = T{r,) (-^ j . (26) 



T 3r 



We take this as an idealized guide for the temperature profile T{t) with proper time. (One 
should note, however, that for temperatures above Tc partonic degrees of freedom contribute 
significantly to the equation of state and thus might modify the here assumed profile sub- 
stantially, if the initial temperature is chosen above the critical temperature.) 

We note that the initial proper time Tq and the dimension D of the expansion are here 
the important parameters determining the dynamics of the expansion: Large D and small 
To lead to a more rapid expansion and cooling. 'Initial' is meant here as the proper time 
To where the partonic gas confines into the mesonic degrees of freedom and before the roll- 
down. We thus choose the critical temperature Tc as the initial temperature. (We will also 
later comment briefly for cases where we have chosen higher initial temperatures.) For the 
(unknown) initial volume V{tq) we will take Vq = 10 — 200 fm? as a reasonable range, which 
implies a spherical initial domain of radius r = 1.4 — 3.6 fm. (Later we will see that varying 
the initial volume will not lead to a major change in the final results within our model.) 

In order to make a statistical analysis we need to sample the initial configurations $a(To) 
and $a(To) at the initial temperature in a systematic manner. As demonstrated in the 
last subsection we let the chiral fields propagate at thermal equilibrium for sufficiently long 
hypothetical times at the initial temperature in order to generate a consistent ensemble of 
possible initial configurations for $„ and $a. The main assumption here is thus then the 
hypothesis of (nearly) perfect thermal equilibrium for the initial chiral order fields before 
the possible roll-down period. 

In [0 the average and statistical properties of individual solutions of the above Langevin 
equations ( [25| ) within the markovian approximation (cf. ( pOD and (|2T|)) have been studied 
with the emphasis on such periods of the time evolution when the transverse mass /x_l 
becomes imaginary and therefore an exponential growth of unstable fiuctuations in the 
collective fields might be expected. It was found that for different realistic initial volumes 
individual events lead to sometimes significant growth of fiuctuations. For the quantification 
of the resulting strength of the coherent pionic zero mode fields and as an experimentally 
more direct and relevant quantity we consider in the following the effective pion number 
content of these chiral pion fields. In the semi-classical approximation this number is given 
by 

n. = \n. [^\r) + ^^\r))^ Vij) . (27) 

This expression can be most simply obtained by considering the energy density of the zero 
mode e^;fc=o = 1/2 (m^Tr^ + vf ). As vr^(T) will be proportional to 1/V(r) at the late stage of 
the evolution after the roll-down period, nT^ij^ then becomes constant at late proper times 
when the effective pion mass [i^ relaxes towards its physical vacuum value. This constant 
number will be extracted from the simulations as the total pion number freed from the DCC 
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decay. For this effective pion number one crucial point is then how large the evolving volume 
V{t) of the DCC domain has increased when the pion oscillations have emerged. 

We will now first employ our model to understand the effect of the dissipation and the 
possible role of memory effects on the evolution. We then further investigate within different 
scenarios the statistical distribution in the resulting pion number (pTj) and will propose a 
new signature of stochastic DCC formation based on the cumulant expansion. 

At this point one might indeed worry why we only consider the k = zero mode and 
not also some other long wavelength pionic excitations, which should also experience some 
unusual amplification according to the general wisdom of DCC formation. From a principle 
point of view our model could be worked out or generalized to take into account also some 
more long wavelength modes. The cutoff momentum should then be taken as < fee ~ 
a/A/^ ~ 400 MeV to account for the pionic modes who could possibly become unstable and 
thus amplified. From the power spectrum shown in the work of Rajagopal and Wilczek 
one notices that even within the drastic quenched situation of instantaneous cooling only the 
lowest discretized momentum mode becomes dominantly amplified, whereas the next higher 
lying pionic modes only show some moderate behaviour. In the more physical situation the 
inclusion of a thermally generated mass term A/2T^(r) in the effective potential will cut 
down even further the low momentum range for possible unstable modes, i.e. kc <^ vA/vr. 
Furthermore also the volume V{to) of an initial domain as chosen by us (at T ^ Tc) is much 
smaller than in P], with a radius between 1.4 - 3.6 fm. Hence, in such a quantized picture of 
a finite volume only a few Fourier modes except the zero mode could really become unstable. 
We therefore expect that only the pionic zero mode can predominantly be amplified. 

III. DISSIPATION: MARKOVIAN VS NONMARKOVIAN 

DESCRIPTION 

In this following section we address on a quantitative level the possible differences be- 
tween the full nonmarkovian treatment and the markovian ('instantaneous') approximation 
for the dissipative ( ^0|) and fluctuating dynamics (|2T|) within the Langevin model. 

The exact nonmarkovian functional r[$] of ([20| ) at a given temperature T and plasmon 
mass rup has been worked out in appendix ^ (As also stated in the appendix we only 
consider in the present investigation the contribution of thermal scattering to the dissipation 
functional, i.e. the part denoted as 71 in the appendix.) As ellaborated in [Q and stated 
in the equation (|T^ the appropriate markovian limit for a sufficiently weakly dissipatively 
interacting system results in the on-shell dissipation or viscosity coefficient rj = r{uj = rrip), 
i.e. (^). For the nonmarkovian dissipational functional we therefore consistently choose 
for the plasmon mass rUp the temperature dependent transversal mass fj,± = mT,(T) of the 
right upper picture of Fig. |. Besides of evaluating a history dependent memory functional 
r[$] to treat the full nonmarkovian dissipative dynamics, as a further complication one 
also has to face the problem of how to simulate colored (i.e. non-white) gaussian noise for 
the fluctuating forces in order to be consistent with the underlying fluctuation-dissipation 
relation ([T3|) or (|2T|) . Our strategy for achieving a numerical realization of colored gaussian 
noise is briefly outlined in appendix ^. With this we can then numerically solve the full 
nonmarkovian equations of motion. The Langevin equations ([T7|) or (|25D are then solved 
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for both cases by means of a standard third order multistep scheme, the Adams-Bashforth 



method [40 



In a strong couphng theory hke the hnear cr-model and also for instable situations en- 



countered in describing possible DCCs the magnitude of the soft modes y |0P might vary 
sufficiently fast so that no dominant oscillatory frequency of the fields does occur and thus 
the Markov approximation should not hold. This gave the motivation for this particular 
study. As it turns out, and as we will argue in the following, however, for situations (and 
thus appropriately chosen parameters for D and Tq), where large and experimentally signifi- 
cant DCC can occur, the distinction between the two cases becomes more or less irrelevant. 
One can then incorporate the numerically much simpler markovian treatment. On the other 
hand, to the best of our knowledge, our study represents the first numerical treatment of 
nonmarkovian Langevin equations in thermal quantum field theory and might certainly be 
of importance for other related topics, e.g. in the description of phase transitions in cosmo- 
logical settings by means of Langevin equations |^ . 

A general expectation for the possible difference is that the rate of thermalization, i.e. 
how fast the considered relevant modes do approach their thermal equilibrium properties 
within the heat bath, might be substantially affected. This is best and most straightfor- 
wardly demonstrated for very simple classical examples like Brownian motion of a diffusive 
particle or oscillator. For a more systematic investigation in this respect we refer to a future 
publication |^, where also the difference between 'weak' and 'strong' dissipative Langevin 



behaviour for diffusive processes will be discussed. 

Referring to our present model it is certainly interesting to study how fast the order 
fields can move (or 'diffuse') towards their equilibrium properties discussed in the previous 
sections [11 B| and [11 C| . (A somewhat similar study for simple markovian dissipation has been 
previously carried out in |^.) In Fig. |^ we show for various cases the relaxation of the 
ensemble averaged a field {{a)){t), being initially distorted by hand, towards its equilibrium 
value {a)eq{T) (compare Figs. ||and|^) in a surrounding heat bath at fixed temperature. As 
constant volume we haven taken Vq = 100 fm?'. In the two upper figures we haven chosen as 
initial values $a(t = 0) = (0,7ri Ri (|$|)eg(T), 0, 0) and <i>a(t = 0) = (0,0,0,0), i.e. an initial 
distortion of the chiral zero mode fields in one particular pion direction along the effective 
finite temperature dependent chiral circle a^ + vr^ = (|<l'p)eg(r). In the upper figure the 
situation is depicted at the critical temperature Tc-, whereas for the middle figure we have 
taken T = 80 MeV . For this investigation we consider 10^ independent simulations for taking 
the ensemble average. For both cases the averaged a field follows a damped oscillation along 
the effective chiral circle. For the markovian simulation one sees that the relaxation towards 
equilibrium goes in accordance with the value of the dissipation coefficient (^) depicted in 
Fig. 1^. The nonmarkovian evolution now shows a slightly less damped relaxation towards 
the equilibrium value, the difference being more pronounced for the lower temperature. 
Qualitatively one can understand this behaviour by comparing the frequency spectrum of 
the dissipation kernel r(a;) (its reduced form is shown in Fig. |T9|) with the on-shell damping 
coefficient used in the markovian approximation T{uj) = T{mp) = rj. This spectrum has its 
maximum in frequency more or less exactly at the on-shell frequency, so that simulation 
carried out within the markovian approximation will result in an effectively larger damping 
and thus faster relaxation, since the effectively contributing frequency modes $(u;) of the 
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motion $(t) in the full treatment are less damped (for uj ^ rup) than those in the markovian 
approximation. 

Another interesting example is shown in the lowest part of Fig. |^. Here we consider 
the relaxation of the order parameter, initially being distorted to its vacuum value, towards 
its equilibrium value {a){Tc) = (7f)(Tc) = at the chiral phase transition without exphcit 
symmetry breaking. Here both the effective masses fi± and fi\\ of the chiral fields vanish, 
so that the effective potential does not posess any quadratic term. Looking again on Fig. 
p!9| one would expect from the behaviour T{uj — »• 0) ^ of the dissipation kernel for low 
frequencies that within the nonmarkovian treatment the relaxation towards equilibrium will 
be much prolonged. This trend can certainly be seen from inspecting the figure. However, 
the nonlinear effective 0^ potential drives the initial relaxation comparable to the simple 
markovian treatment. A significant and steadily increasing reduction of the relaxation rate 
sets in only at later stages of the evolution, when the effective potential really becomes flat. 
The complete relaxation within the nonmarkovian scheme shows thus a highly nonlinear 
behaviour. 

We now go over to discuss the possible differences for the dynamics of the order parameter 
including the simple D-dimensional expansion and cooling scenario discussed in section [HP 



in view of possible DCC formation. As the characteristic quantity we concentrate on the 
effective final pion number n^ of (pTf). Potential DCC pionic modes are driven by the initial 
as well as the intermediate fluctuations experienced in the evolution. 

Generally it is clear that dissipation will subsequently diminish potential large pionic 
fluctuations and thus also decreases the strength, i.e. the pion number, of the potential 
DCC candidate. Only a sufficiently fast expansion and cooling, where the expansion and 
cooling rate is comparable or larger then the experienced damping rate, can counterbalance 
the effect of dissipation on the heat bath. To start to be more quantitative let us consider 
first the markovian description. One has to compare the Raleigh damping term D/r (the 
effective 'Hubble' parameter) with the dissipation or viscosity coefficient rj. Both associated 
terms in the equations of motion (pSj) will diminish the amplitude of any pionic fluctuations 
being buildt up during the roll-down period. On the other hand the effect of the Raleigh 
damping on the pion number content n^ is exactly counterbalanced by the volume dilution 
(p^). n^ being buildt up during the roll-down period can thus physically only be decreased 
by the 'true' dissipation experienced from the heat bath. Whether this dissipation can act 
substantially depends on whether the damping coefficient is comparable in magnitude to the 
Hubble parameter 

D _ D /T\3/^ 
r To \TcJ 

In Fig. ^ we compare f]{T) with the Raleigh coefficient D/t(T) for 3 set of parameters 
of dimensionality D of the expansion and initial proper time Tq. This serves as a rough 
illustration how fast the expansion has actually to be for any potential DCC candidates to 
appear. For some reasonable choices of D and Tq one can see that the Raleigh damping D/t 
will be sufficiently larger than r], at least for later temperatures below about 70 MeV. For a 
sufficient fast expansion D/r will be much larger than t], so that the dissipation due to the 
interaction with the heat bath can not have any tremendous effect on the potential DCC 



candidates except for a slight hindrance on the evolution. The important thing during the 
roll-down is that the fluctuation due to the noise will be large and can eventually enable 
a large disorientation of the order parameter. For moderate or slower expansion, however, 
when both damping coefficients becomes comparable in magnitude after the roll-down period 
even for later times, the dissipation will lead to an additional strong reduction for the pionic 
fluctuations and thus for the pion number, making DCC formation physically impossible. 

In order to support these qualitative arguments we calculate the average pion number 
(i.e. the sum of the pion number of each individual event divided by the total number of 
events) and the pion number of the 'most prominent' event within 10^ independent events by 
solving the markovian Langevin equation ( p5D and compare those with the result obtained 
by solving the same equation but without the damping term rj and the fluctuating noise. 
(The thermally distributed initial conflgurations are the same for both cases.) The most 
prominent event is meant here and in the following sections as the one where the flnal pion 
number is the largest within the generated, flnite ensemble. The 'most prominent' event 
is at flrst, of course, of no direct statistical signiflcance. The error of its occurrence for a 
flnite ensemble will indeed be very large. We explicitely show it for the reason to simply 
see what maximum magnitude in the pion number is possible within a flnite total number 
of generated events within one particular chosen ensemble. 

The calculations are performed for different parameters D and tq. Table | shows the 
results. For a discussion and possible motivation for the various parameters and their actual 
physical relevance we refer at this point to the next section |I^. Here we want to stress 
that the results of table | conflrm our arguments: For the relative slower expansion the 
dissipation due to the interaction with the heat bath destroys any possible large pionic 
oscillations and therefore leads only to a small total pion yield. In contrast, for a fast 
expansion the dissipation has only a minor influence on the DCC formation. For these cases 
the damping coefficient t] is indeed rather small compared to the Raleigh coefficient D/t. 

Now we can answer our primary question: Is there any difference between the full non- 
markovian treatment of the dissipative dynamics compared to the approximate markovian 
treatment on the possible formation of DCC. From our findings at the beginning of this 
section we expect that the effective damping experienced by the memory effects within the 
nonmarkovian case is moderately, but not significantly diminished. ('Memory' indicates 
that the earlier stages of the evolution influence the present motion of the order parameter.) 
The answer is 'frustrating' and simple. For a moderate expansion the pion yield from the 
DCC decay will increase compared to the markovian treatment, but only slightly. In any 
case for such a situation the possible pion yield obtained within the simulations are too 
small to have any experimentally relevant consequence! On the other hand, for a sufficiently 
fast expansion (, which might be speculative or not to be realized in a relativistic heavy 
ion collision,) and for which more prominent DCC candidates will show up (compare the 
next section), the memory effects of the treatment of the dissipation and noise are not of 
particular significance for the final pion number distribution. 

As one particular example we show the outcome of a simulation for both cases in Fig. 
p. We take D = 1 and tq = 0.5 fm/c to simulate a somewhat moderate expansion. ( The di- 
mensionality parameter D = 1 for longitudinal expansion simulates a rather slow expansion. 
On the other hand the here chosen value of tq is very small so that the initial cooling and 
expansion after the onset of the phase transition is still rather fast. This value is definitely 
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too small to be realized in nature. Typically one expects a few fm/c for the onset of the 
phase transition. In the next section we will see that only a. D = 3- dimensional expansion 
can lead to any prominent DCC candidates for reasonable choices of Tq. Therefore one 
should not consider this present example as a physically relevant scenario. Its purpose is 
merely to be an example which does indicate some differences between the nonmarkovian 
and markovian treatment.) In the upper part the sampled distribution of the final pion 
number is shown within 10^ independent events. In the lower part the individual trajec- 
tories (cr, 7ri)(r), i=l,2,3, for the most prominent event out of each ensemble are depicted. 
Clearly these trajectories represent the ones expected intuitively for a true DCC event. This 
intuitive picture is further strengthened when examing Fig. [1^, where the evolution in time 
of the transversal mass ^i±{t) and the effective pion number nT^{T) are depicted for the most 
prominent candidate within the markovian simulation. One clearly recognizes that for this 
candidate the evolution starts with an unstable situation where the pion fields and thus also 
the pion number will be amplified significantly in the very first stage. The subsequent minor 
decrease in the pion number is then attributed due to the further experienced damping. 

The point to make here is that the pion yield of the most prominent event as well as 
the average pion number in the ensemble are somewhat larger within the nonmarkovian 
simulation. However, in sake of the variety of choices for the parameters D and tq and 
the associated wildly differing outcomes (compare e.g. tables |I| and fT[), a modification as 
presented here is only of minor importance. 

We certainly have now to answer what this torturing enterpise for achieving a simu- 
lation of nonmarkovian Langevin equations was good for! In general dissipation as well 
as the associated noisy fluctuations are non-local phenomena in time leading to a memory 
functional over the past history of the system for describing the dissipation as well as to 
a finite correlation in time of the noise. The more phenomenological markovian and white 
noise approximation are generally used in one way or the other as their numerical realiza- 
tion are considerably more simple. Typically such an approximation is justified in a loose 
sense when there exist a clear separation of timescales between the slow degrees of freedom 
under consideration and the ones integrated out. At first sight this is not really given for 
our situation, though our investigation shows that one might indeed work with the much 
simpler markovian approximation. It is also easy to imagine that such a separation is not 
given either for a variety of interesting problems in other areas of physics where one wants 
to describe the effective dynamics of a system in terms of only a few 'relevant' degrees of 
freedom. We therefore believe that our investigation and in particular the numerical real- 
ization of nonmarkovian Langevin equations with colored noise is of general and principle 
relevance for similar problems of classical or quantum dissipative systems in other areas. 
In addition, our extensive discusssion here underlines the importance of understanding the 
certainly complex dissipative nature of the chiral phase transition in more detail. We believe 
that our 'choice' for describing the dissipation for the pionic fluctuations of the zero mode 
with the surrounding heat bath is motivated by an intuitive physical picture. If, on the other 
hand, one can show that the experienced dissipation for the pionic (transversal) fluctuations 
is in fact much stronger, then there is no chance at all for any DCCs to be formed in heavy 
ion collisions! 
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IV. STOCHASTIC FORMATION OF DCC 

Although the inclusion of dissipation as discussed in the last section ^ gradually de- 
stroys on general grounds any possible large oscillations of the coherent pionic fields during 
and after the roll-down period, there still should be a chance for a particular large pion yield 
originating from some 'appropriate' initial fluctuations of the order parameter and as well as 
from the subsequent fluctuations experienced. If the initial fluctuations are large, and if the 
subsequent expansion during the roll-down of the order parameter is sufficiently fast, there 
should indeed be a considerable probability for a long-time large disorientation of the chiral 
fields away from the a-direction during and shortly after the roll-down period. This would 
then lead to a particular large final pion yield. As emphasized already in |jT6|, individual 
statistical events will lead to sometimes significant growth of pionic fiuctuations. In this 
sense the formation of a particular 'large' DCC, i.e. with a sizeable amount of low mo- 
mentum pions being emitted, can follow some unusual distribution to occur because of the 
special stochastical and nonlinear dynamical nature with a possible, temporarily onsetting 
instability. To anwser this question of how often particular events might occur with some 
unusual large pion yield, we investigate in the following the distribution of the pion number 
for different DCC scenarios which differ from each other in the cooling or/and the sampling 
of the initial fluctuations. We will see that the distribution in the flnal pion number takes 
a nontrivial and nonpoissonian form, at least for the more speculative scenarios or param- 
eters employed where one might expect larger DCCs to occur. By means of the cumulant 
expansion of the resulting distributions we will then show that the higher order factorial 
cumulants are even still moderately large when allowing for an additional and incoherent 
realistic background of low momentum pions. Therefore these unusual fluctuations might 
indeed be observed experimentally and thus provide a very interesting new signature for a 
nonequilibrium chiral phase transition and the associated formation of DCC. 

As a crude estimate for the maximum soft pion number to occur from the decay of a DCC 
one can think of a 'true' DCC where the chiral order flelds 'circle' around with the maximum 
amplitude as (vr^ H — ^vf ) < /^ along the chiral circle a"^ + n"^ ~ /^ at some intermediate 
stage f after the roll down in the evolution (compare with the lower part of Fig. |^). (Due 
to the ongoing expansion in our model the amplitude will then subsequently decrease due 
to the experienced Raleigh damping, so that at late times the chiral fields will then only 
fiuctuate around the vacuum value {a) = f-j^.) This will result in a coherent pion number 
density of n-„/V ~ /^m^/2 ^ 0.08 /m~^. For the total pion number the crucial question is 
then how large has the evolving volume V{f) of the DCC domain increased when the pion 
oscillations have emerged. 

At this point we should give also another rough estimate of how many low momentum 
pions should emerge out of decaying DCC in order for a chance of experimental detection. 
In a relativistic heavy ion collision at RHIC one typically expects around 1000 pions being 
produced per unit in rapidity. On this 'background' one has to look for a peculiar and unusual 
enhancement in the pion spectrum at low transverse momentum to identify possible DCC 
formation. It is thus clear that the number of emitted pions out of a domain should be 
somehow comparable to the number of background pions for a particular small window of 
low transverse momentum. The expectation is that one should have a surplus of at least 50 
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pions stemming from a DCC per unit rapidity in a window oipt < 200 MeV in order to allow 
for a promising detection ||4^ (see also the schematic Fig. |16|). This number should thus 
serve in the following as a rough guide. We will come back to the experimental detection 
possibilities at the end of this section. For our above estimate for the maximum number 
of pions out of a domain this would mean that the intermediate volume V{f) has to have 
increased up to a value of about 10'^ fm^ when the order parameter has reached the chiral 
circle. This again implies to consider (or demand) a rapid expansion, i.e. to consider (D=)3- 
dimensional expansion and sufficiently small initial time tq, as already demonstrated in the 
last section (see table |) and which first was emphasized in other studies [0,^. 



A. Different scenarios 

In the following we will present numerical results for the formation of DCC, i.e. the coher- 
ent amplification of the pionic chiral order fields resulting in a final pion number (^) being 
effectively emitted by the domain, for various parameter sets and also for four somewhat 
different scenarios (see also |]TB[)- 

The first scenario we want to discuss is the 'normal' one already described in section 

EB 

Langevin or 'annealing' scenario [|T^,|T^]: Like the dynamical calculations in the last section 



m the initial configurations of the chiral fields (and their velocities) are sampled statistically 
for an assumed thermal equilibrium at the initial temperature T{tq) = T^ and an initial 
volume V{tq), thus covering a (nearly) complete set of possible initial thermal conditions. 
For r > To the subsequent evolution of the chiral order fields for each individual realization 
of the sample is then described by the (markovian) equations of motion (^) within a D- 
dimensional scaling expansion according to ( ^61) and (|2^) . 

Particular examples are listed in the tables | and y and in the Figures |^-|T2| for various 



parameters D, tq and Vijo). As outlined at the beginning of section [II D| one expects that 
the chiral phase might set in at proper times tq ~ 3 — 7/m/c. Inspecting table | one 
recognizes that employing a,D = loiD = 2 dimensional scaling ansatz either the average 
outcome ((ri7r)) or also the outcome for the most prominent candidate of the sample are 
unacceptable small for experimental detection. Only for the D = 3 case (see the tables 
and the Figures listed above) individual and unusual events might occur for a small initial 
proper time tq < ?>fm/c and which might be detectable. This is the situation for a very 
rapid expansion and cooling as noted the first time by Randrup |Tl|. We note, however. 



that the average ((^vr)) is still only moderate even for this rapid scenarios, i.e. ((^vr)) ~ 3 — 4, 
and thus also unacceptable small (see e.g. the upper part of Fig. |T2|). As already stressed in 
]TB| , for an experimental identification this would imply to look for (rather) rare and unusual 
strong fiuctuations on an event by event analysis in certain rapidity and low pt— windows. 
On the other hand one can clearly recognize from the outcome that the original annealing 
picture proposed by Gavin and Miiller [0 and assuming there a rather moderate expansion 



and cooling (tq ~ 7 — 10 fm/c) does not work as the final pion number is by far too small 
(confer table 0). Experimentally significant DCCs cannot happen for this picture according 
to our calculations. 
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Inspecting table || more closely, one recognizes the at first sight maybe somewhat paradox 
behaviour, that the average as well as the pion yield for the most prominent candidate of 
each numerically generated ensemble do not show any strong sensitivity on the chosen initial 
volume V{tq). This result one can at least qualitatively understand as follows: The initial 
fluctuations of the chiral fields depend on the initial volume as discussed in section [11 C| (see 
also Fig. 1^). For a smaller initial volume the initial fluctuations become stronger. Hence 
there is a larger probability for the order parameter to start to evolve against the positive 
a-direction into the 'backward hemisphere'. (For this see e.g. the two examples shown in 
Fig. ^ All 'more prominent' candidates do show a time evolution for the chiral fields akin 
to the ones depicted there.) For such a case the order parameter has to turn back during 
the roll-down so that period of large pionic fluctuations would be prolonged. This then gives 
raise to a larger disorientation of the order parameter. On the other hand, however, the 
volume of a DCC domain at the freeze-out time is accordingly smaller for the initial volume 
being smaller. Referring to ( P7D both trends seem to nearly exactly counterbalance each 
other and hence leading to this peculiar behaviour. 

We now turn to present a few results obtained for three other, more speculative scenarios: 

• quench scenario: The initialisation at T = Tc follows completely analogous to the Langevin 
scenario. However, when switching on to the evolution for t > tq including the volume di- 
lution (|2^), we demand that during the expansion the term AT^/2 of the effective potential 
in the equations of motion (^) is being omitted as well as also the dissipative term and 
the noise in order to mimic an abrupt occurence of the zero temperature vacuum potential. 
Consequently the dissipation and the fluctuation vanish during the roll-down and the oscil- 
lation of the order parameter. Within this scenario we try to simulate somewhat the picture 
proposed in @,0], where it is assumed that the effective potential below T^ changes quasi 
abruptly to the vacuum potential for T = 0. However, the initial conditions, on the other 
hand, are sampled at thermal equilibrium at critical temperature Tc. We believe that this 
picture represents a strong idealization and probably is not likely to happen in an ultrarela- 
tivistic heavy ion collision. Due to the abrupt cooling there is likely more instablity for the 
order parameter allowing for stronger final fluctuations. 

• modified Langevin scenario: This scenario differs from the Langevin scenario only in the 
sampling of the initial configurations. For the sampling we neglect the explicit chiral sym- 
metry breaking (i.e. H = 0; see Fig. P). On the other hand, for the evolution at r > tq 
we employ the same equations of motion (j2^) including the explicit symmetry breaking. In 
the modified Langevin scenario the initial fluctuations are stronger than for the other two 
scenarios since the most probable initial value of the order parameter is centered around 
(((a)) = 0, ((tt)) = 0) and the effective potential for this case is more flat. Hence the possibil- 
ity for the order parameter to start its evolution towards the backward hemisphere is more 
likely to occur. One might argue that the use of the initial conditions prepared within this 
picture is inconsistent within the linear sigma model with a physical pion mass. When dis- 
cussing Fig. H we noted that the phase transition resembles a smooth crossover. From QCD 
lattice calculations, however, one knows that the chiral transitions happens much sharper 
within a very narrow window close to T = Tc. This means that it might very well be that 
the order parameter will (strongly) fluctuate around zero near the critical temperature as 
mimiced by the present realization of the initial conditions. With this in mind one might 
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consider this present scenario even more realistic than the Langevin scenario following the 
simple minded linear a-model. 

• modified quench scenario: The initial configurations of the order parameter are sampled 
as in the modified Langevin scenario. The dynamical evolution corresponds to the quench 
scenario. 

In Fig. [Tl| we depict the distribution P{n.,^) of produced pions logarithmically within 



10^ events within the four different DCC scenarios. As parameters we choose D = 3 and 
To = 7 fm/c [|T3| and V{tq) = 100 /m^, i.e. still only a rather moderate expansion. As 



expected, the pion yields in the modified scenarios are larger than in the normal scenarios for 
the most prominent event as well as for the average. A comparison of annealing and quench 
scenarios, both with finite and vanishing pion mass (for generating the initial conditions) 
reveals that the most productive DCC events would lead for this set of parameters to a 
few (6-8 in annealing scenario with finite pion mass), to a moderate number (20 - 40 in 
annealing scenario with zero-centered initial conditions or quench with massive pions) or 
to about 140 long wavelength pions (in quench scenario with initial conditions generated 
by massless pions), respectively. The final results, of course, majorly depend on how fast 
the effective cooling and expansion proceeds, i.e. on the value of the initial time Tq and 



thus the overall initial Hubble constant D/tq (see also Fig. |T^). In general one finds that 
for sufficiently fast expansion individual unusual strong fiuctuations of the order of 50 - 
200 pions might occur in all the four scenarios, although the average number ((n^)) of the 
emerging long wavelength pions only posesses a rather moderate (and likely undetectable) 
value of 5 -20. 

For a direct comparison we depict a poissonian distribution 

P(n) = ^e-% (28) 

n! 

where the mean value n is equal the averaged pion number ((rijr)) obtained numerically for 
each sample. With the chosen parameters of Fig. ^ the distribution of the pion number 
for the Langevin scenario is indeed still similar to a simple poissonian distribution. (As 
mentioned above, for such a slow expansion the coherent pions produced from a DCC decay 
would be washed out by the background of incoherent pions and thus could not provide any 
signature.) However, for the other three cases the final distribution does noHoUow an usual 
poissonian distribution. This represents a very important outcome of our previous [0 and 
the present, more detailed investigation! Fluctuations with a large number of produced pions 
are still likely with some small but finite probability! In principle, an ensemble averaged 
description of potential DCC formation carried out within the mean field approximation, 
as presented in the various literature, can not account for such fiuctuations and thus has to 
fail at some point. We remark further that also the so called isospin ratio signal is close to 
that expected for a DCC event. 

To demonstrate this interesting behaviour of strongly nonpoissonian fiuctuations even 
more pronounced, we show in Fig. the pion number distribution obtained within the 
Langevin and modified Langevin scenario for a rather fast expansion {D = 3 and tq = 
3 fm/c). This parameters are in line with the ones used in other studies 0,|11|,^,Q. Both 



distributions differ strongly from their corresponding poissonian distributions. The averaged 
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pion number ((72,7)) are 3.9 and 18.5, respectively, and are both comparable to the values 



obtained within the quench and modified quench scenario of Fig. |TI|. The appearance of 
particular events with very large pion number (more than 200) is hereby attributed to the 
initial fiuctuations and the ones experienced during the roll-down periode. 

This special and unusual statistical distribution (obtained within the modified Langevin 
scenario) will be further analyzed in the next section |IVB| . We want to note at this point, 



that Bjorken and collaborators ^^ had speculated that the final distribution in the pion 
number will arguably go beyond a gaussian (or poissonian) distribution even when the DCC 
fiuctuations are generated by a gaussian distribution for the initial condition parameters. 
This 'expectation' is what we have demonstrated now. The initial conditions (compare 
Figs. I and H) follow, more or less, a gaussian distribution, whereas the final occuring 
distribution in e.g. the pion number strongly deviates from a gaussian behaviour for the 
(assumed) nonequilibrium situations, where the occuring DCC phenomena could be consid- 
ered as experimentally be detectable. We also like to mention that Krzywicki and Serreau 
had recently found in a somewhat similar setting, following the model of |^8[, that the 



so called enhancement factor for the final fiuctuations also will follow some unusual and 



nonpoissonian distribution |^ (see also |1T6|). 

There also had been the conjecture in the literature that the long wavelength amplifi- 
cation of the pionic fiuctuations is not really driven by the 'true' DCC phenomenon, but 
actually could be attributed to a parametric resonance behaviour driven by the late and final 
oscillations of the coherent cr-field [^,^. This alternative idea we can at least qualitatively 



address. In Fig. |1^ we show the statistical distribution of the final number for the 'cr'-quanta 



by means of an analogous expression as (pTj), i.e. 

na = ]:m, f (a(r) - ((t),,,)^ + -^a^r) 1 V{t) , (29) 




relaxing to a constant value in the late oscillations of the cr-field in longitudinal direction 
around its vacuum value. The scenario and parameters chosen are the ones for the very 
pronounced situation of the lower part of Fig. 0. On the average about ((?^o-)) ~ 3 cr- 
particles are produced, but within the sample also some events with more than 30 a-particles 
can occur. Due to the potential vacuum decay a — > tttt (with a width on the order of a few 
hundred MeV) the existence of these cr-quanta would result on the average in 6 additional 
pions or, for the more pronounced events, up to more than 60 additional pions. This is a quite 
reasonable number, however, it is still considerably smaller than the direct pions stemming of 



the true DCC as seen from the lower part of Fig. |I2|. So, if these late a-oscillations really do 
exist, because of energy conservation, the amount of pions being produced out of them, either 
thinking in a pertubative way as a result of an individual decay of a cr-quantum or within 
the nonpertubative mechanism of parametric resonance, the numbers of produced pions is 
found to be significantly less than the direct ones of 'true' DCC in all simulations carried 
out. We therefore are tempted to conclude that parametric behaviour is not as efficient 
compared to the 'true' DCC phenomenon. On the other hand, we want to stress here, that 
these late cr-oscillations obtained in our simulations are actually a caveat of our model. As 
already pointed out at the end of section |IIB| , exactly because of the possible decay mode 



a -^ TTTT, one has in principle to account for an additional temperature dependent dissipative 
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term in longitudinal direction for the evolving order parameter. An inclusion would in fact 
then accordingly continously decrease these 'radial' oscillations because of the decay into 
pions [p3| , p4| . 



As a last investigation we consider the possibility that one might prepare the initial 
conditions for the fluctuating chiral fields at some higher initial temperature Ti{Ti) ^ Tc 
within the Langevin scenario. The order parameters are then centered more or less around 
zero []TT|. Switching on to a rapid 3-dimensional scaling expansion one intuitively would 



expect that the chiral fields still would fluctuate around zero when the system cools down at 
and below the critical temperature Tc and thus providing somehow similar initial conditions 
like in the modified Langevin scenario. Accordingly one would expect a more dramatic 
yield in the pion numbers comparable to the one obtained within the modified Langevin 
scenario. It turns out that this is not the case. We find that the final yield follows more 
closely the result of the Langevin scenario with the initial conditions sampled at T = T^, if 
the parameters chosen for the initial time Tj and initial volume K(tj) are adjusted in such 
a way that they exactly coincide with tq and V{tq) for the standard Langevin case when 
T{t = To) becomes Tc. In Fig. |1^ we show the statistical distribution of the a-field and 
its temporal gradient at r = Tc = 3 fm/c for the case when the system was prepared at an 
initial temperature Tj = 300 MeV. The final pion number distribution (not shown) looks 



more or less identical to the upper part of Fig. 12, i.e. to the corresponding Langevin 



scenario being prepared at T = Tc. From Fig. |14| one notices that indeed the a-field is still 
more centered around zero than within the standard case (compare with Fig. ^. However, 
the distribution of the gradient a has shifted towards a nonvanishing positive value because 
of the forward drift experienced by the explicit symmetry breaking term. This shift in the 
later distribution at T Ri Tc then effectively influences the outcome especially in the final 
pion yield distribution in a counteracting way compared to the naive expectation. 

Summarizing this subsection let us highlight once more the main outcome of our in- 
vestigations: If in a heavy ion collision it will come to the formation of an experimentally 
detectable DCC domain - say with at least 50 low momentum pions - it very likely has to be 
a rare event with the average yield {{n.,r)), stemming from the zero mode fluctuations from a 
single domain, still being considerably smaller than 50. A dedicated event-by-event analysis 
is then unalterable. If, on the other hand, nature is more 'obliging', it might also be that 
{{nn)) ~ 50 (as in some very speculative quench scenarios), there exist again with some finite 
probabilities on the percent level some events which contain a multiple in the number of 
pions compared to the average. One can consider those particular events as really unusual 
'pion bursts'. Also for such a situation a dedicated event-by-event analysis is definitely de- 
sirable. For both cases, a possible detection of unusual fluctuations would provide nontrivial 
evidence for the formation of DCCs and the existence of the chiral phase transition. Of 
course, all this is speculation, as also the whole issue of possible DCC formation is. Which 
of the scenarios or assumed parameters are realized in nature one does not know. A slow or 
moderate expansion of the system within the Langevin scenario, which one may consider as 
the most physical one, will indeed not result to any verifiable signal. In the next section we 
will discuss in more detail on the statistical nature of the unusual distributions found and 
on their experimental detection possibilities. 
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B. Critical dynamical fluctuations 

In the last subsection |1V A| we have demonstrated that the probabihty distribution in 
the number of the coherently emitted soft pions from the DCC decay is nontrivial and non- 
poissonian for a sufficiently fast expansion in the various scenarios presented. The statistical 
facets of these unusual probability distributions are what we want to explore in more detail 
in this last subsection. 

Although in fact the distributions in the pion number from a DCC might be realized 
as such, one very likely can not prove these directly from the experimental measurement of 
the unusual pion number abundances, as there are much more pions emitted independently 
at the late stage of an ultrarelativistic heavy ion collision [|^. Since one expects that the 
emission of the soft pions would be affected most significantly compared to the moderate or 
high (transverse) momentum pions, one has to consider to allow for a low momentum pt-cut 
in the data to enhance significantly (and sufficiently) the signal to background ratio. With 
this at hand we will then show in the later part of this subsection that indeed the unusual 
fiuctuations might still be clearly visible and thus provide a very interesting and new event 
by event signature for DCC formation to be analyzed via a cumulant expansion in the (to 
be) measured low momentum pion number distribution in a given rapidity interval. 

In order to account for the true higher order correlations of the statistical distribution 
Pijij^) we consider as a characteristic tool an expansion in factorial cumulants 9m- For a 
rather brief introduction and some further properties and analysis we refer to appendix 
0. The factorial cumulant 9m of order m = 1,2,3,... represent the non-trivial statistical 
m-point correlations of the distribution. 

In [J^l it was stated that so called bin- averaged factorial reduced cumulants for higher 



than two (i.e. m > 3) are consistent with zero when analyzing the particle multiplicity of 
(lighter) nucleus-nucleus collisions (and contrary to hadronic collisions). From this fact Elze 
and Sarcevic then motivated to describe the occurring multiparticle density fiucutations in 
such reactions by means of a (gaussian) three dimensional statistical free field theory p7|| . 



and suggested the conservative view that no first or second order phase transition should be 
imphed as long as there is no compelling evidence in the data. Our situation, of course, is 
different as we (have to) assume a rapid chiral phase transition to occur in order to mimic 
the formation of DCCs. 

With the probability distribution of the soft pion number obtained numerically within 
our model we can calculate the factorial moments and subsequently the factorial cumulants. 
In appendix |C| we have stated explicitly the first six factorial cumulants expressed via the 



usual factorial moments of the distribution. In Fig. |I5] we show the first six reduced factorial 
cumulants 9m/ < n^^ >™ in a logarithmic scale for different expansions simulated by D = 3 
and varying tq within the modified Langevin scenario (see also the lower part of Fig. |12|) . 



Here < n.,^ >= ((^tt)) = 9i denotes the average pion number of the corresponding distri- 
bution. Each distribution, except one, was sampled by 10^ independent events employing 
the markovian equations of motion. For the one remaining distribution (with tq = 3 fm/c) 
10^ events were generated in order to estimate the possible error. One immediately realizes 
the striking behaviour that the higher order and reduced factorial cumulants 9m/ < n-,^ >™' 
with m > 3 are clearly nonvanishing and in fact show an exponentially increasing tendency, 
as for each distribution the higher order reduced cumulants lie more or less on a straight 
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line in this logarithmic representation! Comparing the results for the two distributions for 
To = 3 fm/c, where one was generated with a sample of 10^ events and the second one with 
a sample of 10^ events to account for higher statistics, we can estimate the error for the 
higher lying factorial cumulants for a sample of 10^ events to be still in the order of factor of 
two. This can easily be traced back to the obvious fact that the higher factorial cumulants 
depend most sensitively on the tail of the numerically generated distribution with large mul- 
tiplicity n. On the other hand the general trend of exponentially increasing reduced factorial 
cumulants is not affected by the higher statistics [[ This behaviour suggests a special sort 
of dynamical scaling as the (at least) higher order factorial cumulants approximately take 
the form 

e^ ^ ae"" < n >'" , (30) 

where a and a denote constant numbers, depending on the parameters chosen for sampling 
of the distribution, and < n >= {{n.^)) just represents the average number of pions of the 
distribution P{n.,t). With this asymptotic form (PD| ) for the factorial cumulants 9m (now 
assumed also for m = 1 and m = 2) one can actually invert the expansion and find the 
corresponding distribution P{n) giving raise to such characteristic factorial cumulants. This 
is briefly worked out in appendix y. The resulting distribution, reflecting for the higher 
order factorial cumulants, is given by a shifted Poisson distribution 

n' 

Pin) = T^e-'^ (31) 



in' 



'V 



where n = {< n > e") ■ n' and n' = 0, 1, 2, . . .. As generally a is positive and a small number 
- typically from the above slopes one has a -C 1 - and ((n7r))e" is some multiple of the 
average number {{n.„)), the deduced distribution (|3lD provides a nice intuitive and intriguing 
picture for the unusual events: Occasionally, a semi-classical 'pion burst' with pion number 
n^ = (((?T-7r))e°') ■ n' is being emitted for some special events. These represent rare events 
as the distribution in n' follows a standard Poisson distribution sharply peaked at n' = 0. 
Such rare and unusual events are then in fact quite similar to the Centauro candidates 0! 
We do not want to push this interpretation too far, as smaller deviations from the straight 
exponential fit and, of course, the two lowest factorial cumulants are not considered. Yet 
we believe that this interpretation provides the right intuitive way of describing the unusual 
strong fluctuations in the tail of the distribution. 

At this stage one might indeed ask for the physical origin of such a peculiar and scaling- 
like behaviour of the fluctuations. Here we can provide at present no definite answer, as 
we can only rely on our numerical findings. For a given ensemble of initial configurations, 
the stochastic approach presented in this work results in an ensemble of widely differing 



^ On our suggestion, an expansion in higher order cumulants in the distribution of the so called 
enhancement factor ^O; as given in [^], shows exactly the same tendency of exponentially increas- 
ing reduced cumulants |48]. As the model there is different from ours, this repeated finding points 



towards some 'universal' behaviour. 
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solutions. Since a gaussian initial distribution in the fields under the time evolution of a 
quadratic Hamiltonian always stays gaussian, we believe that the unusual final fluctuations 
in the present case originate due to the particular nonlinear evolution. In principle, the 
occurrence of some sort of peculiar scaling behaviour in higher order factorial moments, is 
known (or speculated) for quite a time to show up in the multiplicity fluctuations stemming 
from a quark-hadron phase transition (in hadron-hadron or heavy ion collisions) described 
within a simple phenomenological Ginzbug-Landau framework [^. In this respect our 
findings underline the necessity to learn more about the possible onset of a phase transition 
by a careful study of final multiplicity fluctuations. 

On the other hand one also clearly recognizes that the second order factorial cumulant 
increases drastically compared to any 'usual gaussian' second order cumulant on the order 
of the first order cumulant and thus defines a much broader distribution. This increase in 62 
is due to the fact that many trajectories of the sample enter temporarily the instable region 
with ij,\ < 0. In fact for the more dramatic cases we have 62 ^ Oi. In appendix |C| we briefly 
show that for a situation, where 62 > Oi, there exists no simple statistical distribution which 
can be expressed solely in terms of the first two factorial cumulants. This thus signals again 
that the statistical nature of the distribution is highly non-trivial! 

It remains to shed light on the possibility whether such unusual fluctuations can indeed 
also be reflected in the cumulant expansion based on the data measured in real ultrarela- 
tivistic heavy ion collision experiments. In the real world one expects a huge 'background' of 
pions not coming from the decay of a 'large' DCC domain, as already outlined at the begin- 
ning of this section |^. To illustrate such a background we show in Fig. |l^ a schematic and 
qualitative expectation of a single event of the transverse momentum spectrum of charged 
pions within some definite rapidity interval including a single hypothetical and sufficiently 
prominent DCC candidate. Such a spectrum has been schematically redrawn from a single 



simulated event of background pions to be expected at RHIC energies |^3|. A 'large' DCC 
domain would eventually enhance the number of soft pions in the pion spectrum at suffi- 
ciently low momenta (see Fig. |16D. The authors of [^] provide a detailed analysis that if 



allowing for a low momentum p^-cut of pt < 200 MeV in some small and definite interval of 
rapidity (of order one) the expectation is that one should have a surplus of at least 50 pions 
stemming from a DCC per unit rapidity in such a window for a possible 'direct' observation. 
Even in such a small window, however, if supposedly large DCC domain occurs, there will 
be still a background of 'normal' pions of the order of 50 in average. Therefore the inherent 
fluctuations of the background pion number in low pt makes it rather difficult to find out 
a clear trace of the DCC formation in the soft pion enhancement within one event. One 
then has to go to an appropriate statitical analysis for discovering possible unusual fluctu- 
ations. More importantly, as we have stated in the last subsection, 'larger' DCC domains 
are more likely to be some rare events. We thus want to pursue in the following by means 
of the cumulant expansion whether there is a possibility to look experimentally for unusual 
fluctuations when allowing for some additional incoherent and simple fluctuating poissonian 
source producing also low momentum pions and thus providing the background. 

By now there have been two experimental investigations to look for DCC events either in 
heavy ion reactions at the CERN-SPS [^] or in p — p-collisions at the Tevatron at Fermilab 
|5T| . Both programs had so far a negative outcome in their searches. This might still well be 



due to the fact that up to now no analysis employing a sufficiently small low momentum cut 
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has been carried out. In addition, also a wavelet-type analysis, as originally been proposed 



by Huang and coworkers |£2|, might further help to look for the occurrence of unusual events 
or - in respect to our present work - of unusual fluctuations in sufficiently small rapidity and 
momentum windows. There exist also other clever suggestions how to filter for the DCC 



events, see e.g. 53 



At this step we want to provide a rough estimate in what range the typical rapidity 
interval is to be expected for the pions to be emitted out of a single DCC domain. The 
spherical (D=3-)scaling expansion ansatz in proper time was chosen to mimic for the rapid 
expansion. In strict terms such a scenario is fueled by everlasting sources and thus should 
break down at some later decoupling time as the whole collision of two heavy ions does 
last only a finite time. Before "freeze-out" the domains are separated from the outside or 
exteriour vacuum by the surrounding and expanding matter. This deficiency of everlasting 
sources can be circumvented by a mapping of the idealized 3-dimiensional boost-invariant 
evolution to quasifree, truncated sources evolving in normal time at some decoupling time as 



shown by Bjorken and coworkers [44]. Such a truncation of the evolution modifies somewhat 



the final momentum spectrum of the emitted pions ^^. In any case, within the idealized 
scenario evolving solely in proper time at least at the beginning of the evolution, a simple 
estimate for the rapidity interval is given by 

1, /l + t^c 



A,.ilnfi±^) 
' 2 \l-VcJ 



where Vc = , , ^° = and r(ro) = i^Vir^Yl'^ . For the parameters employed (tq = 2 — 

7 fm/c and V{tq) = 10 — 200 fm?) this estimate implies a rapidity interval of Ar^ ^ 0.2 — 1 
for the low momentum pions to be emitted, in agreement with general expectation. 

Suppose now that we have the following situation: The soft pions are coming from either 
a DCC domain or, independently, from an incoherent ('chaotic') source (background). Fur- 
thermore we assume the emission of the incoherent soft pions follows a standard poissonian 
distribution with the mean value < n,r >p, i-e. 

As the resulting (factorial) cumulants in the independently combined pion number distri- 
bution are additive (see appendix P), the reduced cumulants can thus simply be written 
as 






m 



< n^ >™ (< ra^ >^ + < n^ >^nc^rn 

where 'c' denotes the coherent emission by a DCC state and 'inc' the incoherent emission 
by the background source. For m > 2 the cumulants related to the incoherent pion source 
do vanish by assumption of a Poisson distribution. (We note that if there are more than 
one single domain contributing within a considered rapidity and momentum window, and 
if these are truely uncorrelated, the respective cumulants of each independent source would 
then again simply add up for the combined pion distribution.) Fig. |1^ depicts the resulting 
reduced factorial cumulants ( PBj) obtained for a single domain simulated with a fast expansion 
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{D = 3, tq = 3 fin/c, modified Langevin scenario, compare with lower part of Fig. ^ and 



Fig. 15) and superimposed by tlie inclusion of a background source with different mean 
values < tZtt >p ranging from 20 to 200 incoherent additional pions. The last numbers can 
either be seen simply as basic uncertainty and/or also as a result of lowering the pt-cut. 
The additional poissonian source basically lowers all the reduced factorial cumulants with 
rra > 1: As in its form (p2D , P^p'^{n) has no factorial cumulants 9m with m > 1, the combined 
reduced factorial cumulants become smaller accordingly. However, the higher order ones for 
m > 3 are still appreciably large if the background mean pion number is less than about 
100, especially with increasing number m. If we consider as an example the situation that 
< Utt >p is about 70, it shows that for a slow expansion tq = 7 fm/c the reduced factorial 
cumulants of higher order m > 3 are still very small (< 10^^). This basically reflects 
the suppression of the few coherent emitted pions compared to the large background. In 



contrast, however, for a fast expansion tq = 3 fm/c, as depicted in Fig. ^ where 'large' 
DCC states are more likely to occur, the reduced factorial cumulants of higher order are in 
the range 1 ~ 10 and thus should be clearly visible and detectable. 

We thus find that for sufficient fast expansion the reduced higher order cumulants are 
still in the order 1 ~ 10, although the number of incoherently emitted pions might in 
average be 3 ~ 4 times larger as the averaged number of DCC pions. (In our particular 
last example we have {{nT,))Dcc = 18.5 compared to < n^ >p= 70.) We therefore conclude 
that an experimental analysis by means of the higher order factorial cumulants for the 
low momentum pion number distribution provides a well-suited indication for the possible 
existence (and to some lesser extent also for the identification) of any DCC formation on an 
event-by-event analysis. Event-by-event type analysis for getting additional new insight in 
the underlying physics of heavy ion collisions (e.g. the process of thermalization) has become 
quite popular over the last two years (see e.g. [Q and references therein). In this respect 
our work can be considered as a special and ultimate scenario of what to expect in case of 
a rapidly ongoing chiral phase transition associated with possible DCC state formation. 

V. SUMMARY 

In the present work we have elaborated in detail within an idealized, but microscopi- 
cally motivated semi-classical Langevin description on the statistical facets of the formation 
of possible disoriented chiral condensates during and after the onset of the chiral phase 
transition expected to occur in ultrarelativistic heavy ion collisions. Within the Langevin 
treatment of the standard linear cr-model, one can simulate, on an event by event anal- 
ysis, the possible evolution of various DCC scenarios in a rather transparent form. Our 
main focus and objective has been to understand the physical role of dissipation and noisy 
fiuctuations on the DCC phenomenon. The advantage of the presented approach is that 
in contrast to common mean-field treatments, which can only bring about a deterministic 
description for the (ensemble) averaged evolution, it allows for any possible branching of the 
dynamical trajectories being especially important in the instability region. Our Langevin 
picture is based on microscopic input, although one can interpret the presented approach 
more intuitively also in the spirit of the phenomenological Landau- Ginzburg description of 
phase transitions. Our ideas could also be taken over for situations advocating a first order 
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transition within the hnear cr-model [^, in order to study for such parametrizations of the 



effective temperature dependent potential the influence of dissipation and fluctuation on the 
evolution of the order parameter inside the nucleating and growing bubbles. 

The model, originally being first proposed in |jl6|, is based on the very assumption that 
the high-momentum particles ('hard' fluctuations) of the chiral fields constitute a heat bath 
which behaves locally thermalized in the expanding system. The interaction of the nonequi- 
librated 'soft' chiral fields with this surrounding heat bath then gives raise to their stochastic 
and semi-classical evolution of Langevin type. Our main conception is that the order param- 
eter as well as the pionic fluctuations before and after the onset of the chiral phase transition 
still interacts (dissipatively) with its surrounding of thermal (or 'hard') pions, which then 
results in large and tremendeously differing fluctuations during the evolution. Furthermore 
we have concentrated solely on the effective dynamics of the collective zero mode (order 
parameter and pionic fluctuations). We have argued, that, if at all, the zero mode pionic 
fluctuations become most unstable during the roll-down period and thus are the ones being 
predominantly amplified for realistic initially small sized and separated expanding domains. 
The overall picture of possible DCC evolution resembles the one proposed by Bjorken and 
coworkers P,^. 

As a first application we considered the finite size fluctuations of the order parameter and 
the chiral pionic fields for a given volume and temperature, resulting in a further smoothening 
of the crossover behaviour around the critical temperature within the employed linear cr- 
model. The Langevin description provides a powerful and simple tool for generating in a 
systematic and efficient manner a canonic sample of statistically possible configurations at 
a given temperature T. As a reasonable (or minimal) assumption we consider the order 
parameter and the chiral fields to be likely thermally distributed when the phase transition 
during the later expansion starts to occur. This sampling of possible initial configurations 
contrasts to the ad hoc guesses for the initial conditions made in many of the previous 
works on DCC physics and thus enables us to investigate characteristic statistical properties 
of DCC formation. 

We have then concentrated on the dynamical evolution of one single domain during and 
after the onset of spontaneous chiral symmetry breaking at the later stages of the initially 
very hot system expected to occur in ultrarelativistic heavy ion collisions. Because of the 
collective expansion at these later stages, the temperature will subsequently drop below 
the critical one, and smaller, originally chirally restored domains (assumed to be indepen- 
dent being separated spatially and in rapidity) start to form together with a thermahzed 
background of (quasi-)pions and possibly other hadronic excitations within the respective 
expanding subsystem. A D-dimensional scaling expansion was employed to account for the 
collective expansion resulting in an additional Rayleigh or Hubble like damping term within 
the stochastic equations of motion. 

We stressed the important issue of the physical effect of dissipation on the pionic fluctua- 
tions for any possible DCC evolution. For the quantification of the resulting strength of the 
coherent pionic zero-mode field and as an experimentally more direct and relevant quantity 
we considered the effective pion number content n^^ (via eq. (^Tj)) of the emerging final os- 
cillations in the chiral pionic fields. The dissipation kernel has been calculated by means of 
a standard finite temperature field theory technique and is directly associated to the inverse 
thermal scattering rate of the soft mode on the thermal particles. Our analysis clearly shows 
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that the (rapid) expansion, i.e. the Hubble damping term, has to be at least as efficient 
in order to compensate for the true dissipation. For a larger dissipation coefficient r] the 
final yield in the pion number would be correspondingly smaller, as the dissipation damps 
accordingly faster any large DCC like pionic fiuctuations which have possibly emerged after 
the roll-down. Although our estimate for the dissipation close to the critical point is inspired 
by physical arguments, further understanding of the certainly complex dissipative nature of 
the chiral phase transition is crucial: If one can show that the experienced dissipation for 
the pionic (transversal) modes close to the transition point is in fact much stronger than the 
one we have employed, then there is definitely no chance at all for any DCC signals to be 
seen in heavy ion collisions. On the other hand, as emphasized at the end of section [11 B| , in 
the deeply broken phase much below the critical temperature the associated dissipation for 
the pionic fiuctuations will certainly be reduced by the additional chiral cr-meson exchange 
for the 71 — TT-scattering amplitude as compared to our employed estimate. For the deeply 
broken phase then also the decay a ^ vrvr becomes possible giving potentially raise to a 
much stronger dissipation for the fiuctuations in longitudinal direction. All this might very 
well at least quantitatively effect some of our results and conclusions presented concerning 
the possible survival of DCC states and would require an even more involved and detailed 
calculation. 

In addition, we also have described (in appendix ^) how to numerically realize colored 
noise in order to treat the underlying dissipative and nonmarkovian stochastic equations 
of motion. In general dissipation as well as the associated noisy fluctuations are non-local 
phenomena in time. This, to the best of our knowledge, is the flrst numerical treatment of 
nonmarkovian Langevin equations in thermal quantum fleld theory and might certainly be 
of relevance for other related topics. 

In the last section we have then given a comprehensive numerical study for the possible 
formation of DCC, i.e. the coherent ampliflcation of the pionic chiral flelds, for various 
parameter sets and also for four somewhat different scenarios. It shows, as pointed out 
the flrst time by Randrup [0, that a rather rapid expansion is mandatory to have any 



signiflcant chance for obtaining 'large' DCCs which then might lead to some experimental 
consequences. On the other hand, our analysis has provided the at flrst sight more pessimistic 
view, that even then, a DCC event has very likely to be an unusual and rare event. The 
average characteristic ((^tt)), i.e. the average number of low momentum pions being emitted 
of the flnal pionic modes, shows only a moderate behaviour, which then should result, on 
the average, in a mild increase of the transversal low momentum spectrum in the pions. As 
we have argued, such a mild increase is probably tremendously difficult to observe directly 
and unambiguously from the average momentum spectrum of pions. 

However, the statistical distribution P{nT^) of emitted pions shows a striking nonpoisso- 
nian and nontrivial behaviour. There exist within some still flnite probability some rare and 
unusual events which contain a multiple in the number of pions compared to the average. 
As pointed out in the subsequent analysis of the statistical nature of such distributions, one 
should indeed interpret those particular events as unusual and semi-classical 'pion bursts' 
similar to the mystique Centauro candidates . This result suggests a very important con- 
clusion: If DCCs are being produced, an experimental flnding will be a rare event following 
a strikingly, nontrivial and nonpoissonian distribution. A dedicated event-by-event analysis 
for the experimental programs (e.g. the STAR TPC at RHIC) is then unalterable. 
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We clearly have to say once more, that, of course, all the above conclusions represent 
speculation, as also the whole issue of possible DCC formation is. Which of the scenarios 
or assumed parameters considered in our work are realized in nature one does not know. 
Any slow or moderate expansion of the system will indeed not result to any verifiable signal! 
In addition, because of the various approximations and idealized scenarios considered, our 
work should not be seen as directly comparable to any experimental data. In this respect, 
our final last theoretical conjecture, which we now want to summarize, has to be seen as a 
fascinating and experimentally possible guideline for a future analysis of the pion spectra 
to be taken at RHIC or already taken at CERN-SPS, if DCC like phenomena occur in 
ultrarelativistic heavy ion reactions. 

For any meaningful experimental identification our results imply to look for rare and 
unusual strong fluctuations on an event by event analysis in certain rapidity and sufficiently 
low {pt) windows. The further analysis of the unusual distribution in the pion number asso- 
ciated to a rapid chiral phase transition we have invoked by means of the factorial cumulants 
dm-i which represent a powerful tool, well-known in the analysis of final multiparticle fiuc- 
tuations in high energy hadronic reactions. We have found the striking behaviour that the 
higher order and reduced factorial cumulants 9m/ < ^vr >™ with tti > 3 show an abnormal, 
exponentially increasing tendency. This we consider as the most important outcome of our 
extensive investigation. In addition, we also found that the second order factorial cumulant 
62 increases dramatically compared to any 'usual' gaussian distribution, thus characterizing 
a much broader distribution. This broadening refiects the fact that many trajectories of the 
sample have entered temporarily the instable region. In addition, we have allowed that on 
top of the pions emerging from the decay of collective pionic modes a further incoherent and 
poissonian background source of low momentum pions might in fact overshadow or even 
completely wash out these striking characteristics. As it turned out, however, the reduced 
higher order factorial cumulants are still of the order 1 ~ 10, if the number of incoherently 
emitted pions is already in average 3 ~ 4 times larger than the average number of DCC 
pions. 

We therefore strongly advocate that an analysis by means of the higher order cumulants 
serves as a new and powerful signature to identify any unusualities associated with potential 
DCC formation. Of course we are aware that our last analysis assumes, that within each 
window in momentum and rapidity, where the experimental analysis is considered, a DCC 
like phenomena with conditional probability equals to one has occurred. This might not be 
the true case. However, we believe that our suggestion for future experimental analysises 
is in fact rather 'simple' to carry out and represents most likely the only way to find (any) 
evidence for unusualities in the low momentum pion spectra. If such an analysis turns out 
to be negative, there is probably no other chance to look for the DCC phenomenon. 
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APPENDIX A: THE DISSIPATION KERNEL r(K = 0,a;) 

In this appendix we evaluate the frequency dependence of the dissipation kernel r(k = 



0,uj) used in section |T| for studying the nonmarkovian dissipative evolution of the chiral 
fields. 

For this we had employed the sunset diagram from standard $^-theory generalized to the 
present 0(4)-case, which will result in a different numerical coefficient and will be specified 
at the end of the appendix. 

In the $^-theory the respective dissipation kernel for soft modes (with |k| < kc) inter- 
acting with the hard modes (with |k| > kc) separated by a momentum cutoff kc is given 
as 

r(,..)^!MJ^f. --^"^"---) |. |k|<.„ (Ai) 

2uj \ u 



where the memory kernel A4 introduced in ||23| reads 



TT 4 f <fqid?q2 1 



24 Jkc (27r)° co'ico'2'^3 

{[(l + ni)(l + r22)(l +^23) - nin2n^]5{u - ui - c^2 - ^3) 
+[ (1 + ^1)722(1 + "n-s) - '"■1(1 + n2)n^]5{ijj - ui + uj2 ~ cus) 
+[ni{l + n2)(l + ns) - (1 + ni)n2n3\6{uj + lui - iU2 - UJ3) 
+[^1^2(1 +213) - (1 +ni){l + n2)ns]6{uj + uJi + UJ2 - W3) 

+[ (1 + ni) (1 + 722)723 - 7ii7l2(l + n3)]5{uj - UJi- UJ2 + UJ3) 
+[(1 +7ii)7l2723 -72i(l +7i2)(l + n:i)]6{uj - iOi + UJ2 + UJ3) 

+[ni{l + n2)n3 - (1 + 711)712(1 + n3)]6{iu + cui - iU2 + UJ3) 
+[nin2n3 - (1 + 7ii)(l + n2)(l +n3)]6{uj + uJi + UJ2 + ^0-3) }, 

and q3 := k — qi — q2, tUj = cjq., rii = n{uji) = ^^Jt_^ , i = 1,2,3. The dispersion 



relation for the hard modes is taken as ujq = Jq^ + mi, where rup denotes the dynamical 
mass. iM. (k, uj) represents the net absorption rate for soft modes due to the interaction 
vertex of a soft mode with three hard particles. The first and the last term of ( |A^ ) correspond 
to the decay of one soft mode into three hard modes and the inverse process. The other six 
terms correspond to the scattering process s + h ^^ h + h. 

For our study of stochastic DCC formation we constructed an effective model for the 
chiral zero mode fields, so that for the present purpose we take kc = and thus need only to 
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calculate iA4(k = 0,uj). Evaluating the (5-function we will reduce the six dimensional integral 
of (|A2|) to a 1-dimensional integral which we then treat further numerically. In principle this 



task had already been performed by Wang and Heinz |^6[ investigating the 2-loop resummed 
propagator for hot $^-theory. However, repeating the steps in their tedious derivation we 
found out that some particular kinematic boundaries of the integration variables were not 
extracted correctly. In the following we sketch the main strategy and then state the final 
result for the considered dissipation kernel. 

It is easy to see that iA^(k, uj) is antisymmetric in uj and therefore r(k, uj) is symmetric. 
We thus only have to consider the case for a; > 0. In this case the contributions from the 
4th, 6th, 7th and the last line of eq. ( [A2| ) are identical to zero. In addition one sees that the 
contributions from the 2nd, third and the 5th line are the same. Moreover one can convince 
oneself that the yield of the respective absorption processes just gives a common factor er 
compared to the respective emission processes due to the standard detailed balance relation 
for systems at thermal equilibrium. With these observations we have 

tMiO,u) = ^j^ (e^ - l) fd'qid'q2 -^— x 
24 (27r)° ^ ^ J UJ1UJ2UJ3 

{3(1 + 721)722^3 6{uj + uji- UJ2- UJ3) 

+nin2ns 6{uj - lui - uj2 - ^3)} ■ (A3) 

We now outline our strategy by manipulating the first integral in (|A3|) which corresponds 
to the emission rate of the scattering process h + h -^ s + h. This integral can be reduced 
to a 3-dimensional integral 

g,{uy.=3 fd\d% ^^ + '''^'''''^ S{u + 00, -002- 00,) (A4) 

J UJ1UJ2UJ3 

= 24iT / dqidq2 dt q^q^ d[uj + uji - 002 - OJ^,) , 

Jo J-l UJ1UJ2UJ3 

where t = cosO and 9 denotes the angle between qi and q2. For the energy conservation 
stated by the ^-function one has 



UJ + UJ1-UJ2 = CJ3 = Jqf + gl + 2gig2^ + rnl (A5) 



due to the momentum conservation qi + q2 = qs- ( [A5| ) represents the kinematical constraint 
among the variables gi, q2 and t. In order to determine from this equation for a given 
frequency u the kinematic boundaries for gi, q2 and t we take its square yielding 



iF{q2, t) :=) ^''_^''f /' - -^q2 t = ^1+^ = u;2{q2) ■ (A6) 

Taking qi as the most outer integration variable we now consider it as a fixed constant 
and concentrate first on the variables ^2 and t. The left side of (|A6|), which we define as a 
function -F(g2, t), represents a straight line in ^2 with different inclination for different values 
of t. All straight lines for different t G [—1, 1] cut at g2 = 0. Then the solutions of (|A6|) 
for a fixed gi (and given a;) are the points where the bundle (in t) of straight lines F{q2,t) 
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cuts 072 (g2)- There are three cases to distinguish and which are classified by the position of 
F{q2,t) at g2 = 0: (case I) F{q2 = 0,t) > rrip] (case II) < F{q2 = 0,t) < rrip; and (case III) 
F(g2 = 0,t) < 0. Fig. illustrates the different situations for the three cases and shows 
the kinematic boundaries of g2 and t. 

We have to remark that the solutions of ( |X6| ) do not necessarily fulfill the original con- 
straint ([A5|). Therefore one has to insert back the solutions into ( |A5| ) and check whether 
they indeed fulfill (|A^). One finds out that the solutions of case III do not obey ( [X5| ) (as 

Ul + Uli — UJ2 < Q) . 

The kinematic boundaries for the integration variable q2 are 

(l2=2\yB{uJ^-qi\ , ql' = - (JB{uj,qi) + qi 

1 / r—, r \ .,' 1 



g2'=2 [-yBiuJ, qi) + qij , 1^' = 2 (^V ^(^' ^1) + ^1 

{uj + uJif \{uj + LUi)^ - g? - 4m2 
where p6| B{uj,qi) = 



(u + uJiY -ql 

In order to fulfill the classification for case I (II) one finds from the definition of F{q2 = 0, t) 
of ( |A6D that the energy u has to be greater (less) than rup. 

One can now get the kinematic boundaries of gi using the fact that the function B{uj, qi) 
should not be negative. For case I one finds that gi has no further constraints. For case II 
gi posesses a lower boundary qf: 



2uj 
Eq. ( |A^ ) can now be stated as 



qT = TTTV (^^ ^ r?7,2)(u;2 - 9^,2 ) 



gi{ijj)=2ATi'^l9{mp — Lo) / dqi I dq2 / dt + 6{uj — mp) / dqi I dq2 / dt 

Ja"' J a^^ J — 1 Jo J a'L^ J — 1 



Ql g2 0{'^ + Ui - UJ2 - UJ'i) ) . (A7) 

(^1(^2^3 J 

By suitable substitutions for the integral variables, 

Q1Q2 ^i 

dt — > duJs = dt , and dqi, dq2 -^ dUi, (if/2 with Ui := e~~ z = 1, 2 , 

one can carry out the integrations over t and g2. The result is 



gi{uj) = 24^T'h{mp-uj)J^ '" dU^ Gi{Ui;q2\qT) 

+e{uj - mp) j^ dUi Gi{Uv, ql\ q^^) \ (A8) 



with 
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and 



Gi{Ui]Si,S2) 
U{s) := exp 



1 



1 



1 - f/i 1 - UiU^ 



■ In 



:i-U{s^)){U{s2)-UiUj 
:i-U{s2)){U{si)-UiU^) 
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T 



U,., :-- 



-ui/T 



The further evaluation of the second term in (|A3|) corresponding to the emission rate of 
the (inverse) off-shell decay process h + h + h ^ s follows in analogous but slightly more 
complicated way to the strategy for the first term considered above. 
We state the final result for i7Vl(0,u;): 







-9{uj — Srrij^ 



U{0) 

t 

U{ql) 



U{ql) 
C 
U{qT) 



dUi G2{Ui-qt\qt)+ r'^'dUi Gi{Ui;q},qt 



(A9) 



with 



G2{Ui; si, S2) ■■-- 
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where 



A{u,qi) 



{uj - cji)2 ({to - uif - ql - Ami) 



{u - uJiY - qf 



(We note at this stage that we have obtained different results for the lower kinematic 
boundaries ^2^ and ^2^ ^is compared to the ones given by Wang and Heinz |5^, which were 
there simply set as zero.) 

The dissipation kernel has then the form 



T(0,uj) 



iM{0,uj) 

2u 



9 T ^ fu^ m^ 
1927r3 VT' T 



(AlO) 



where we have defined a reduced dissipation kernel F (^, ^ j which depends only on uj/T 
and iTip/T. In fig. |1^ we show the reduced dissipation kernel for irip/T = 0.1 and irip/T = 1. 
Here 71 and 72 denote the two different contributions to T from the scattering and decay 
process. One recognizes that 72 in fact diverges for a; — > 00. One can furthermore show 
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that 7i has an assyinptotic behaviour ~ Xjuj for sufficiently large uj. In our present study 
concerning stochastic formation of DCC we have neglected the 72 contribution and have 
thus only described the physical dominant scattering contribution of s + /i <-*> /i + /i. The 
dissipation kernel 71(0;) has its maximum very close to the on-shell frequency uj = rup. Its 
shape with frequency is thus effectively governed by two characteristic and independent 
scales, the temperature T and the plasmon mass rrip. The Fourier transform of the reduced 
dissipation kernel, i.e. r(k = 0,t) = 71 (t), for a temperature T = 120 MeV near the critical 
temperature T^ is plotted in Fig. ^. The mass rUp is extracted as the transversal mass n± 
at that given temperature from the right upper picture of Fig. |^. By means of Fig. ^ one 
can estimate that the correlation in time of the kernel extents to about 5 fm/c. 

On the plasmon mass shell u = nip the damping coefficient r{0, rUp) is obviously given 
solely by the scattering contribution: 

r(0,m,) = J 3 (1-^(0)) / dU, G,iU,-4^=q,,q^,^=0). 
IzSn'^nip Jo 

The above integration can be further simpified to 

(1 - f/(0)) J^ dUi GiiUv, qi, 0) = fsp (1 - e-- 
where fsp{x) is the Spence function defined as 

fsp{x) := - dy 



1 y-l 

In the high temperature limit nip <^ T the dynamical mass nip is evaluated by the tadpole 
diagram as nip = g'^T'^ /24. One thus recovers 



which is twice the plasmon damping rate |5y,56 



Generalizing from 0(1) to 0{N) one has an additional (pre-)factor ^^j^ for the dissipation 
kernel and the dynamical mass. Then for the case of the linear 0(4) a-model one also has 
to substitute g"^ -^ 6A. 

APPENDIX B: NUMERICAL REALISATION OF COLORED NOISE 

Here we outline a new numerical method for simulating nonwhite, i.e. colored gaussian 
noise for an arbitrary (nonnegative and symmetric) noise kernel I{uj). 

For this consider the situation of a 1-dimensional Brownian particle interacting with 
its thermal surrounding. The force acting on the particle can be separated into a mean 
(dissipative) part and a random, stochastic part. This random force shall not depend on the 
state of the particle and represents a fluctuating source given as a particular noise sequence. 

Suppose the noise in a time interval [0,T] is composed of a series of pulses which occur 
randomly ||S^. Each puis can be written as a ■ b{T), where 6(r) has a certain uniform shape 
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and a denotes a random height which is allowed to be either positive or negative. Then a 
particular sequence of noise reads 

n 

i{t) = Y.a,h{t~U), tG[0,T]. (Bl) 



In (|BT|) there are one random number n and two random variables, the height a^ and 



the timing (center time) tj. The random number n is now assumed to obey a Poisson 
distribution with the mean value n = nT, where /i denotes the mean counting rate. The 
random distribution of tj shall occur uniformly. Now one has to specify also the statistical 
distribution p{a) of the height a. In the limit of a large number of sufficiently weak pulses, i.e. 
large n and small mean-square value a^ oi p{a), the noise will then approximately be given 
as a gaussian process due to the central limit theorem. Gaussian noise is solely characterized 
by the first two moments 

mm = 

and Uimit'))) = lit, t') = fia^ f ds b{t - s)b{t' - s) . (B2) 

Jo 

For this limiting case one can freely choose the distribution p{a). Most commonly one 
employs a gaussian distribution 

pia) = , — e 2^ . 



For a large /i the distribution of the puis number n has a sharp maximum centered at n. 
Thus we set a fine time scale At and assume that on each time point one given pulse should 
occur. The number n is then fixed by T/At, and the timings tj of occurence of the pulses are 
also fixed. It remains to find out the form 6(r) which is related to the correlation function 
I{t,t') of (|B2|) . If the correlation function is a ^-function then it is easy to show that 6(r) 
will be as well proportional to a 5-function. For this case one denotes the so constructed 
fiuctuation C,{t) as white noise. In more general case the noise is called colored noise. 

For the simulation of gaussian colored noise we now assume that 6(r) has a symmetric 
shape within some time interval [—A, A]. Outside this interval 6(r) is taken to be zero. The 
noise ^(t) is a stationary process for t G [A,T — A]. It means /(t, t') = I{\t — t'\) = I{t) 
where t = t ~ t'.(For t G [0, A] and t G [T — A,T] there are switching on/off artifacts.) 
Fourier transformation of (p^ ) yields 



Iiu)=fia'\b{iu)\\ (B3) 

From the above equation it follows that one has to demand the Fourier transform /(u;) of 
the correlation function of the noise should be nonnegative. 6(a;) is a real function due to 
the symmetry of &(t). Further we assume b{uj) to be positive. One thus ends with 

b (r) = ^— G(r) for r G [- A, A] (B4) 



with G{t):=I p^)e- 
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For the case of white noise with unit strength one has I{t — t') = S{t — t') and thus 

b{r) = -^6{t) . 
Then the sequence of white noise can be written as 

iw{t) = Y] ai — —6{t - U) 

n 1 

= E ^ci,6{t - U) , (B5) 

4 = 1 V/^ 

where Oj can be sampled according to a gaussian distribution with a unit mean-square value. 
Having fixed the pulse number n by T/At = fiT, the mean counting rate is /i = 1/At. 
Furthermore we approximate the 5-functon as 

*"•' = { t : l7l 
Then the white noise at each time step can be simply generated as 



At 

Coming now back to the construction of a colored noise sequence, we find that ^(t) can 
be obtained by an integral of the history of a particular white noise sequence folded with 
the uniform pulse 6(r) of (|B4|), i.e. 



e(t) =J2a^ h{t - u) 
j=i 

= yaj dt' h{t-t')5{t' -ti)= I dt' h{t-t')yai5{t' -U 
f^i Jo Jo t^i 

= j dt'a^b{t-f)i^{t') 

dt' Git - 1') i^{t') . 



In order to check the reliability of our simulation we calculate numerically the ensemble 
average of the noise correlation and compare it with the given correlation function for which 



we choose the reduced dissipation kernel plotted in Fig. |20|. The result is shown in Fig. |2T 
The depicted average was obtained by 10"^ independently realized noise sequences. 

APPENDIX C: CUMULANT EXPANSION 

We give here a brief reminder of the factorial cumulant expansion for discrete statistical 



distributions [Q . A stochastic number n is fully characterized by its probability distribution 
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P{n), n = 0, 1, 2, . . .. An equivalent and convenient representation of P{n) is given by its 
probability generating function 

oo 

F(l -x) := 5] (1 - x)" P{n) . (CI) 

n=0 

If one defines the factorial moments 0m by 0o = 1 and 

0^ =< n(n- l)---(?2-m + 1) > (m > 1) , (C2) 

then the probability generating function becomes 

oo (_rp\m 

F(i-x) = E ^n^0-- (C3) 

The probability generating function also serves to generate the factorial cumulants 6'^, 
which are defined by 

oo (_rf,\m 

ln(F(l-x)) := Y: ^-^Om. (C4) 

The factorial cumulants 6m represent the nontrivial ('nonirreducible') correlations at order 
m. 

For a Poisson distribution, 

P{n) = -e-\ 

it follows immediately that its factorial moments read 0^ = n™ and consequently we have 
F[l — x) = exp{—xn). Therefore the Poisson distribution is characterized by the vanishing 
of all factorial cumulants except for 6i = n. 

According to the definition (|C4|) the factorial cumulants 6m are generally combinations 
of the factorial moments 0j with {i < m). For our use we list here the expressions up to 
order six, i.e. 



1 

2 = V2 - 01 

3 



6i 
6 

^3 = 03 - 30201 + 20? 

^4 = 04 - 40301 - 302 + 120201 " 60^ 

^5 = 05 - 50401 - 100302 + 20030? + 300^01 - 60020? + 240? 

^6 = 06 - 60501 - 150402 + 30040? " 100^ + 120030201 

-120030? + 300^ - 2700^0? + 3600201 - 1200? . (C5) 

Now we turn to the question of how - or whether it is in general possible - to receive 
the probability distribution P{n), if all the factorial cumulants are given. Using ( |C1| ) one 
obtains 
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Pin) 



d" 



n\ dx'' 
Although one finds out from (|Cl| ) and (|C4D 



-F(l-x) 



(C6) 



n=0 



n) 



F{l-x)\ 



x=0 



the so inverted distribution P{n) does not necessarily be positive for all integers n. Therefore 
the general answer to the question is 'no'. For example assume the simple situation, where all 
6m except 6i and 62 are zero P7| , |5U |. In this case P{n) can serve as a probability distribution 
only for 6*1 > 6*2, i.e. 

n/2 



Pin) 



exp(^2/2 - 01 



n\ 



Hn{i[{e2-e^Y/2e, 



lV2^ 



(C7) 



where Hn{z) is the Hermite polynomial of nth order. For 9i < 62 it is easy to show that -P(l) 
is negative. (This caveat has not been noticed in [^,^.) This means, that for situations, 
where one finds that 62 > 61, there exists no underlying statistical distribution which can 
be expressed solely in terms of the first two factorial cumulants. 

Another interesting example we want to discuss for the purpose of analyzing the findings 
in subsection |1V B| is the case of a shifted Poisson distribution. Here, to some upper limit, the 
resulting factorial cumulants are situated in logarithmic representation on a curve, which can 
be described to a good approximation by a straight line. The shifted poissonian distribution 
we introduce as 



Pin) 



(n 



')\ 



n 



Bn' . 



n 



0,1 



(C8) 



with a and B some positive constant. According to ( PT| ) we have 

aB\ 



F{l-x) 



exp 



B 

E 



-X) 



For small a <^ 1 and large -B 3> 1 the factorial cumulants 9m (with m < B) are approximately 
given by 



On 



a-B' 



(C9) 



The higher factorial cumulants for m > B vanish. 

As a last reminder we consider a combined stochastic process resulting in the discrete 
variable n = ni + n2 by two completely independent stochastic processes PA(rai) and PB(n2). 
The probability distribution PAusin) is then given by 

PAuB{n)= Y. PAini) PBin2) . (CIO) 

ni,n2] 
ni + n2 = n 

With this one finds from the definition ( |C^ ) and ( Pl| ) the factorial cumulants of the joined 
probability distribution PAusin) as 

C^ = o^ + e^. (cii) 

The factorial cumulants of the independently combined variable are additive. 
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TABLES 



(a) D = l 


To (fm/c) 


nTr: the most prominent event/average 


with dissipation 


no dissipation 


1 


7/1.4 


18/3.1 


0.5 


20/3 


34/5 


0.3 


48/6 


60/8 




ih)D = 2 


To (fm/c) 


n-j^: the most prominent event/average 


with dissipation 


no dissipation 


4 


6/1.6 


13/2.5 


2 


22/3 


24/3.7 


1 


75/8.3 


64/8.3 




(c) D = 3 


To (fm/c) 


Tij^: the most prominent event/average 


with dissipation 


no dissipation 


7 


6.5/1.5 


13/2.4 


3 


32/3.9 


28/4 


1.8 


85/9 


65/8.5 



TABLE L The resulting pion yield for the most prominent event and the average, respectively, 
obtained within different expansion scenarios simulated by the special choice of D and tq. The 
calculations have been performed by using the markovian Langevin equation (^). The initial 
volume is chosen as V{tq) = 100 fm^ for all cases. For comparison we neglect the dissipation ('no 
dissipation'), i.e. taking the damping coefficient rj and the noise as zero during the dynamical 
evolution of the order parameter. The results are obtained within an ensemble of 10^ events. 
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I? = 3 and 


t{to) = n 








V^(to) \ TO 


3fm/c 


5fm/c 


7fm/c 


lOfm/c 


10 fm^ 


33.2/2.9 


13.6/1.8 


6.5/1.5 


7.2/1.4 


25 fm^ 


62.8/3.7 


13.4/2.0 


7.0/1.6 


5.3/1.4 


100 fm^ 


30.0/3.8 


13.6/2.1 


6.7/1.6 


6.4/1.4 


200 fm^ 


25.0/3.7 


11.1/2.0 


7.4/1.7 


5.5/1.4 



TABLE II. Pion number of the most prominent event and the average obtained within the 
markovian Langevin scenario for different initial volumes V{tq). The averages are taken over 
2 • 10^ events. The initial proper time tq is also varied to simulate different expansion scenarios. 
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FIGURES 
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FIG. 1. Feynman diagrams contributing to the influence action Sip up to second order 0(5^) 
in S'int- 
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FIG. 2. The temperature dependence of the magnitude of the order parameter |<I>| = \/cj^ + vf^, 
the transversal (pion-like) and longitudinal (cr-like) mass at thermal equilibrium for the physical 
case of a non- vanishing pion mass and the case without explicit chiral symmetry breaking (H = 0) . 
The averages are obtained over an ensemble of 10^ realizations. 
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FIG. 3. The dependence of the friction coefficient r] on the temperature. 



50 



> 

CD 



A 
V 



V=25 fm3 
V=100fm3 
V infinite 




50 100 



150 200 
T [MeV] 



300 




V=25 fm3 
V=100fm3 
V infinite 



250 



300 



150 200 
T [MeV] 

FIG. 4. The temperature dependence of the magnitude of the order parameter ^ at thermal 
equihbrium for different volumes. The averages are obtained over an ensemble of 10^ realizations. 
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FIG. 5. Statistical distribution of the chiral a field and one of the three pion fields for different 
finite volumes. The temperature is taken as the critical temperature Tc. The distributions are 
obtained from an ensemble of 10'^ independent realizations. 
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FIG. 6. Statistical distribution of the chiral a field for different finite volumes for the situation 
without explicit chiral symmetry breaking. The temperature is taken as the critical temperature 
Tr. The distributions are obtained from 10"^ realizations. 
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FIG. 7. Relaxation of the ensemble averaged a field within a heat bath at finite temperature 
for the nonmarkovian and markovian case. The upper and middle part correspond to the situation 
of a physical pion mass, whereas the bottom one corresponds to the case without explicit chiral 
symmetry breaking {m.j^ = 0)- ^^ this case we depict the relaxation of the ensemble averaged 
magnitude of the order parameter. The averages are taken over 10^ realizations. 
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FIG. 8. Comparison of the friction coefficient rj{T) with the Raleigh damping term D/t. 
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FIG. 9. Statistical and normalized distribution P{n.,^) of the emerging final pion number and 
the time evolution of the three pionic trajectories ((cr(T),7rj(T)), i=l,2,3,) for the most prominent 
event within an ensemble of 10^ realizations in both a markovian and a nonmarkovian simulation. 
The trajectories start at the initial proper time tq = 0.5 fm/c. The starting points are marked 
with 'X'. The marks along the trajectories are positioned at time intervals of Ar = 0.21 fm/c. 



56 





300 




250 




200 


-i" 




t/) 


1bO 


U) 




CD 

E 


100 


o 


50 


Q. 





> 




o 


-50 


CI) 




3= 
CD 


-100 




-150 




-200 




- for the most prominent event from 1 

Langevin scenario (markovian) 
D=1, T„=0.5 fm/c, V„=100 fm", T„=T^ 



-10 10 20 30 40 50 60 70 80 90 100 

T [fm/c] 



50 

45 

j_«40 

53 35 

I 30 

c25 
o 
Q. 20 

CD 

> 15 

■I—' 
o 

a 10 

H — 

" 5 





■ for the most prominent event from 1 

Langevin scenario (marl<ovian) 
D=1 , T„=0.5 fm/c, V„=1 00 fm", T„=T^ 



-10 10 20 30 40 50 60 70 80 90 100 

T [fm/c] 

FIG. 10. The time evolution of the effective pion mass /^^(t) and the effective pion number 
n7r(T) for the most prominent event obtained in the markovian simulation of Fig. ^ 
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FIG. 11. Statistical distribution of the final yield in the pion number for four different scenarios 
(see text) within a (D=)3-diniensional scaling expansion. Each simulation has been performed with 
10^ independent events. The initial volume V{to) = 100 /m^ and the initial proper time is taken 
as To = 7 fm/c. The distributions are compared with the corresponding Poisson distributions. The 
averaged pion number in the Langevin, modified Langevin, quench and modified quench scenario 
are 1.66, 2.45, 3.46 and 20.36 respectively. 
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FIG. 12. Statistical distribution of the final yield in low momentum pion number within the 
Langevin (upper figure) and modified Langevin (lower figure) scenario compared with the corre- 
sponding poissonian distribution. A fast expansion is simulated by choosing D = 3 and tq = 3 
fm/c. The upper distribution is calculated within 10^ events, whereas for the lower a sample of 
10^ independent events has been chosen. 
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FIG. 15. The reduced factorial cumulants for m = 1 to 6 for the pion number distribution 
obtained within the (niarkovian) modified Langevin scenario (D=3) for different initial proper 
times. The initial volume is taken as 1^(to) = 100 /m^. The average pion number ((n^)) within 10^ 
events are 54, 18.5, 4.8 and 2.4, respectively, corresponding to tq = 2, 3, 5 and 7 fm/c, respectively. 
In addition the cumulants obtained for a distribution for tq = 3 fm/c within a larger sample of 10^ 
events are also shown to estimate the numerical error. 
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FIG. 16. Schematic and qualitative view of the transverse momentum spectrum of charged 
pions for one single event within some definite rapidity interval including a single hypothetical and 
sufficiently prominent DCC candidate. (This spectrum has been schematically redrawn from a 



simulated event of background pions to be expected at RHIC energies [l43 |. 
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FIG. 17. The reduced factorial cumulants for m = 1 to 6 for the pion number distribution of 
low momentum pions stemming from two independent sources: A distribution stemming from a 
single emerging DCC taken from the lower part of Fig. |l^ (modified Langevin scenario with D = 3 
and To = 3 fm/c) and a poissonian distributed background pion source with different mean values 
{n)p = 20 - 200. 
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FIG. 18. Illustration of the solutions and kineniatical constraints for q^ and t for the three cases 
(see text). For case I there exists one solution for q2 for every t. The kinematical constraints q^ 
and q^ are the points where the lines -F(g2) 1) and F(g2,— 1) cut u}2{fl'2)- For case II there exist 

s' 

two solutions for every t G [tfcn— !]• The kinematical constraints (f^ and q-^ are the cut points 
between ^((72, —1) and W2('72)- For case III all solutions are larger than ^2^. It is easy to show that 
these solutions do not fulfill the equation ([A5|). 
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FIG. 19. The reduced dissipation kernel for nip/T = 0.1 and rUp/T = 1. 71 and 72 are the 
contributions to f from the contributing scattering and decay process (see text). 
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rup = 194 MeV. 
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noise with the given correlation function which is taken as the reduced dissipation kernel of Fig. 
pO| (lower part). The averaging is performed over an ensemble of 10^ noise sequences. In addition 
one exemplaric numerically generated noise sequence dictated by the correlation function is also 
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